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ABSTRACT 

We use very deep UnGTZI multi-field imaging obtained at the Keck telescope to study the evolution 
of the rest-frame 1700A galaxy luminosity function as the Universe doubles its age from Z'-^A to z~2. 
We use exactly the same filters and color-color selection as those used by the Steidel team, but probe 
significantly fainter limits — well below L* . The depth of our imaging allows us to constrain the faint 
end of the luminosity function reaching Mi7oo~— 18.5 at 2;~3 (equivalent to '--^l Mq/jt) accounting 
for both N^^^ uncertainty in the number of galaxies and for cosmic variance. We carefully examine 
many potential sources of systematic bias in our LF measurements before drawing the following 
conclusions. We find that the luminosity function of Lyman Break Galaxies evolves with time and 
that this evolution is differential with luminosity. The result is best constrained between the epochs 
at z~4 and z~3, where we find that the number density of sub-L* galaxies increases with time by 
at least a factor of 2.3 (IIct statistical confidence); while the faint end of the LF evolves, the bright 
end appears to remain virtually unchanged, indicating that there may be differential, luminosity- 
dependent evolution (98.5% statistical probability). Potential systematic biases restrict our ability 
to draw strong conclusions about continued evolution of the luminosity function to lower redshifts, 
2~2.2 and Z'^l.l, but, nevertheless, it appears certain that the number density of ^^--^2.2 galaxies at 
all luminosities we studied, — 22>Mi7oo>— 18, is at least as high as that of their counterparts at ^^^3. 
While it is not yet clear what mechanism underlies the observed evolution, the fact that this evolution 
is differential with luminosity opens up new avenues of improving our understanding of how galaxies 
form and evolve at high redshift. 
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1. INTRODUCTION 

Understanding the formation and evolution of galax- 
ies continues to be one of the most active fields of ob- 
servational cosmology. Over the last decade, advances 
in instrumentation and technique have made it possible 
to select large samples of "normal" star-forming galax- 
ies for direct study at redshifts that correspond to a time 
when the Universe was only a tenth of its present age and 
so study galaxy assembly at a time when galaxies were 
young. Several different approaches for selecting high-z 
galaxies are used, including selection in rest-frame far- 
IR (e.g., Barger et al. 1998; Hughes et al 1998; Blain 
et al. 1999; Eales et al. 2000), near-IR (e.g., Sawicki 
2002), optical (e.g., Thompson et al. 1999; Cimatti et 
al. 2002; Sawicki et al. 2005a) and the UV (e.g., Stei- 
del et al. 1996, 1999, 2003, 2004; Sawicki, Lin, & Yee 
1997; Lowenthal et al. 1997; Giavahsco 2002; Lehnert & 
Bremer 2003; Stanway, Bunker & McMahon 2003; Iwata 
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et al. 2003; Ouchi et al. 2004a). Among these differ- 
ent techniques the Lyman Break Galaxy (LBG; Steidel 
et al. 1996, 2003) surveys have yielded the largest, spec- 
troscopically confirmed samples and the most active and 
detailed follow-up studies. 

Follow-up observations of LBG samples has taught us 
much about the nature of these high-z galaxies, although 
— understandably — such follow-up has so far mainly 
focused on relatively luminous objects at z~3, where 
the samples are largest and spectroscopy is easiest. We 
now know, for example, that LBGs are dominated by 
fairly young episodes of star formation and that they are 
enshrouded by large amounts of interstellar dust (e.g., 
Sawicki & Yee 1998; Ouchi et al. 1999; Shapley et al. 
2001; Papovich et al. 2001; Vijh, Witt, & Gordon 2003); 
that they are associated with massive dark matter ha- 
los (Adelberger et al. 1998; Giavalisco et al. 1998); that 
they have strong, starburst-driven outflows of material 
into the surrounding IGM (Pettini et al. 1998, 2001, 
2002; Adelberger et al. 2003); and that they hkely have 
sub-solar, but not primordial, metallicities (Pettini et al. 
2001, 2002). However, most LBG studies have focused 
primarily on relatively luminous objects (L>L*) at a sin- 
gle epoch (z'-^^S), and comparisons of LBG properties as a 
function of time and luminosity are still in their infancy. 

Studying the properties of galaxies as a function of 
redshift has a straightforward motivation rooted in the 
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fact that observing galaxies at different epochs aUows us 
to study directly their evolution as a function of time. 
The initial search for and study of high-z galaxies has 
been motivated by the desire to find the progenitors of 
present-day galaxies and recent comparisons of galaxy 
populations between different epochs beyond z>l is its 
direct and natural extension (e.g., Steidel et al. 1999; 
Adelberger et al. 2004; Ando et al. 2004; Ferguson et al. 
2004; Papovich et al. 2004). Exploring the time domain 
holds obvious but important attractions. 

The study of high-z galaxies as a function of their lu- 
minosity is less obvious to motivate, especially given the 
fact that high-z studies are observationally expensive 
even for luminous galaxies and far more so for galax- 
ies that are intrinsically faint. Nevertheless, a wealth 
of potential information lies to be discovered in compar- 
isons of galaxies as a function of luminosity. The galaxy 
luminosity function is not a simple power law as could 
be expected from the mass function of dark matter ha- 
los (e.g., Jenkins et al. 2001), but instead reflects the 
imprint of real differences in galaxy formation and evo- 
lution processes. If the luminosities of high-z galaxies 
correlate with the masses of their host dark matter halos 

— as is suggested by some clustering studies (Giavalisco 
& Dickinson 2001; Ouchi et al. 2004b) — then the shape 
of the LF at high rcdshift likely bears the direct imprint 
of star-formation driven feedback as a function of halo 
mass. Additionally, luminosity may also reflect the ef- 
fects of processes such as fluctuating star formation rates 

— whether induced by galaxy-galaxy interactions or by 
other mechanisms — or differences in the properties of 
interstellar dust. In any case, it is unlikely that galaxies 
of different luminosity are just trivially scaled copies of 
each other and so to fully understand the story of galaxy 
formation we must study not just the brightest, obser- 
vationally most accessible members of the population, 
but also their fainter cousins. Such studies are particlu- 
arly attractive because differences in evolution between 
galaxies with different UV luminosities — or, by exten- 
sion, star formation rates — can be expected to point us 
to some of the most relevant mechanisms responsible for 
driving galaxy evolution. Finally, for any reasonable lu- 
minosity function, most galaxies are sub-L*galaxies and 
most of the luminosity in the Universe is contained in 
galaxies below L*. By extension, so is most of the star 
formation and metal production activity. The hitherto 
neglected sub-L* high-redshift galaxies deserve our avid 
attention for all these important reasons. 

One of the most basic descriptors of a galaxy popu- 
lation is its luminosity function (LF). The shape of the 
galaxy luminosity function bears the imprint of galaxy 
formation and evolution processes. The characteristic 
break in the LF seen at both low and high redshift sug- 
gests that galaxies below L* are not simple scaled replicas 
of those above L* but differ from them in more substan- 
tial ways. At present relatively little is known about the 
shape of the faint-end of the LF of galaxies at high red- 
shift, and about the evolution of the high-z LF. What 
studies have been done are limited by small, often single 
fields — such as the HDF — that can be affected by both 
sample and cosmic variance (e.g., Sawicki, Lin, & Yee 
1997; Steidel et al. 1999) or use samples whose fidelity 
has not been well tested with spectroscopy (e.g., Iwata 
et al. 2003; Ouchi ct al. 2004a; Gabasch et al. 2004). 



In this paper we use our large, very faint Keck Deep 
Fields (KDF) galaxy samples to construct the luminosity 
functions of high-z star- forming galaxies over a wide span 
of cosmic time (0.6, 0.8, and 0.8 Gyr from Z'--^4 to z~3 
to z~2.2 to z^l.7, respectively) and reaching to very 
faint hmits {TZ=27, or Mi7oo~— 18, equivalent to star 
formation rate of '--^l Mq/jx). At z~4 and 3, we com- 
bine our samples with LF measurements by Steidel et 
al. (1999) made from shallower, but larger-area surveys 
to study the rest-frame UV-selected galaxy luminosity 
function over up to a factor of a hundred in luminosity. 

This paper is structured as follows. In §|21we briefly 
describe our KDF sample of faint, UV-selected galaxies. 
In §131 we describe the details of how we calculate the LF 
and in § ^ we describe our results and focus on exam- 
ining the possible sources of systematic error that may 
affect them. In § |31 we examine several intriguing evolu- 
tionary trends in the LF. In §|B|we discuss some possible 
interpretations of the observed evolution and also point 
out the potential new approaches to the study of galaxy 
evolution at high redshift that the evolving LF opens to 
us. Finally, in §|7|we summarize our results. As in all 
the papers in the KDF series, we use the AB flux nor- 
malization (Oke, 1974) and adopt r2M=0.3, rL\=Q.7, and 
i?o=70 km s-i Mpc-^ 

2. THE DATA 

Our study of the faint end of the high-redshift luminos- 
ity function uses data from our very deep UnGTZI Keck 
imaging survey — the Keck Deep Fields (KDFs). The 
KDF was specifically designed to explore the evolution of 
the population of very faint (sub-L*) star-forming galax- 
ies at high redshift, including their luminosity function 
and luminosity-dependent clustering. The KDF probes 
the hitherto poorly-explored faint end of the galaxy pop- 
ulation at redshifts z~4 - 1.7 by extending to fainter 
magnitudes the well-known color selection techniques 
used by Steidel et al. (1999, 2003, 2004). The KDF sur- 
vey is described in detail in the companion paper by Saw- 
icki & Thompson (2005; hereafter KDF I), which gives a 
detailed description of the observations, data reduction, 
and selection of star-forming galaxy samples at Z'--^2 - 4. 
Here, we give only a brief overview of its main charac- 
teristics. 

The KDF survey uses the very same UnGTZI filter 
set and color-color selection criteria that are used so 
successfully by Steidel et al. (1999, 2003, 2004) to se- 
lect their high-redshift samples, but reaches a limiting 
magnitude TZiim=27, that is 1.5 magnitudes deeper than 
the surveys by the Steidel group. Because of our use of 
an identical UnGTZI filter set, our KDF data are a di- 
rect and straightforward extension to fainter magnitudes 
of the spectroscopically-tested and well-understood sam- 
ples of Steidel et al. Because of the extensive spectro- 
scopic work of the Steidel team, selection effects, includ- 
ing foreground interloper fractions (which are only a few 
percent), are well known and can be safely used for our 
fainter samples as explained in KDF I. 

The KDF cover a total area of 169 arcmin^ and con- 
sist of five fields (called fields 02A, 03A, 03B, 09A, and 
09B) that were observed separately and so reach slightly 
different depths ranging over TZiim^26.7~27 .3 (50% com- 
pleteness) . The five fields are grouped into three patches 
(patches 02, 03, and 09) that are spatially well-separated 
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on the sky. The division of the KDF into these three 
spatially-independent patches gives us the extremely im- 
portant ability to monitor the magnitude and impact of 
cosmic variance. 

Our ground-based images have typical seeing of 
ensuring that high-z galaxies are unresolved and can be 
treated as point sources. This relatively poor spatial 
resolution is a blessing in disguise as it drastically re- 
duces concerns about galaxy size selection biases that are 
present in observations with better image quality, such 
as HST data. 

Photometry is done in a manner virtually identical to 
that of the Steidel et al. work, namely with object de- 
tection in very deep 7?.-band images and color measure- 
ment through matched 2"-diameter apertures on images 
smoothed to a common seeing. Our UnGTZI filter set — 
the same as that used by the Steidel team — allows us 
to select high-z galaxies in a manner that is identical to 
their brighter samples (see KDF I or Steidel et al. 1999, 
2003, 2004, for details of the selection criteria). The pho- 
tometric completeness of our survey — tested carefully 
using simulations — is similar at Ti=27 to that of the 
Steidel et al. surveys at TZ=25.5. The KDF therefore 
probes a factor of 4 deeper in luminosity than the work 
of the Steidel team. 

To its nominal completeness limit of TZum—^n, the 
KDF contain 427 GT^-I-selected z~4 Lyman Break 
Galaxies (LBGs), 1481 ^7„G7^-selected z-3 LBGs, 2417 
C/„G7?,-selected z~2.2 star-forming galaxies, and 2043 
?7„G7?.-selected z~1.7 star-forming galaxies. 

3. CALCULATION OF THE LUMINOSITY FUNCTION 

We use the effective volume, K//, approach to com- 
pute the luminosity functions. Our approach is virtually 
identical to that used by Steidel et al. (1999) on their 
brighter z~3 and z^A samples. This section of the pa- 
per is devoted to a detailed discussion of the technique 
we use to calculate the LF. We defer the discussion of 
our actual LF results to §^and later. 

The Veff approach to calculating the LF is straight- 
forward. In brief, (1) for each redshift sample (z'~-^4, 
z^3, z~2.2, and Z'--^1.7) we first use simulations to de- 
termine the effective volumes, K//, of the survey as 
a function of apparent magnitude — volumes that ac- 
count for incompleteness due to objects missing from 
the sample. We then (2) combine these with the 
observed galaxy counts to compute the incompleteness- 
corrected number density of galaxies at each redshift as 
a function of apparent magnitude. And finally, (3) we 
convert these apparent-magnitude number densities into 
absolute-magnitude ones to arrive at the LFs. These 
three steps are described in detail in 13.11 and 13.31 
respectively 

We measure the LF at rest-frame 1700 A for two rea- 
sons. Because (1) LF calculation at this rest-frame wave- 
length matches the LF analysis by Steidel et al. (1999) 
using brighter z~3 and z~4 Lyman Break Galaxy sam- 
ples. But also because (2) — as we will see in more de- 
tail in S 13.31 — rest-frame 1700A very closely matches the 
observed-frame /-band at Z'^4, 7?,-band at Z'^3, and G- 
band at z~1.7, thereby nearly eliminating uncertainties 
in fc-corrections. At z^2.2, 1700A is located between the 
G and TZ bands and so for our z~2.2 objects we construct 
composite GTZ magnitudes that more closely match rest- 



frame 1700A than do either G-band or 7?.-band alone. 
Our composite GTZ magnitudes are a simple average of 
the G-band and 7?.-band fluxes, 

G7^ = -2.51og[(10-°-4^ -f 10^o-4^)/2], (1) 

with uncertainties calculated by combining the G-band 
and 7?.-band uncertainties in quadrature. 

3.1. Calculating the effective volumes of the survey 

We must first compute the effective volume of the sur- 
vey, as a function of apparent magnitude, for objects in 
each redshift sample. Because galaxies scatter in and 
out of the color-color selection regions (see Fig. ^ here 
and Figs. 4 and 5 in KDF I) that we use to select our 
high-z galaxy samples, and because this scattering will 
depend on the size of photometric uncertainties, we must 
compute our effective volumes as a function of apparent 
magnitude. Our calculation of effective volumes is ac- 
complished through simulations that we implant into our 
images and then seek to recover artificial objects with col- 
ors and magnitudes representative of star-forming high-z 
galaxies. The description of these simulations is the sub- 
ject of the present section, 13.11 

Modeling the colors of high-z galaxies 

As the first step, we calculate a grid of model colors 
expected of high-z galaxies. We start with model spec- 
tral energy distributions (SEDs) of star-forming galaxies 
from the 1996 version of the Bruzual & Chariot (1993) 
spectral synthesis library. We use the continually star- 
forming models with solar metallicity and Salpeter initial 
mass function (IMF) . We next redden them with a set of 
extinction values in the the Calzetti (1997) starburst dust 
prescription. We then complete the redshift dimension of 
the grid by stretching these SEDs by (1 -I- z) and attenu- 
ating them using the Madau (1995) prescription for con- 
tinuum and line blanketing due to intergalactic hydrogen 
along the line of sight. Finally, we integrate the resul- 
tant reddened, observer-frame model spectra through the 
UnGTZI filter transmission curves to arrive at the pre- 
dicted colors of high-z star-forming galaxies. Some ex- 
amples of model galaxy colors are shown in Fig.^ where 
they are overplotted on top of regions of color-color space 
used to define our galaxy samples. 

The largest influence on the colors of high-z galaxies is 
wielded by, first, attenuation by intergalactic hydrogen 
gas blueward of the Lyman break, and, second, by red- 
dening due to interstellar dust internal to the galaxies. 
Other variables such as age, star formation history, stel- 
lar initial mass function, or metallicity can also play a 
role, but their effects are small in comparison, and espe- 
cially so at the rest-UV wavelenghts that concern us in 
this study. Consequently, we adopt fixed values for most 
of these parameters and explore only how our results vary 
with the adopted reddening and starburst age. 

For practical reasons, we must restrict our choices to 
a limited set of these parameters. We are guided in our 
choice of dust attenuation and starburst age by the ob- 
served values of these quantities in z^3 LBGs. Early on, 
Sawicki & Yee (1998) studied the rest-frame UV through 
optical broadband photometry of 17 spectroscopically- 
confirmed z^3 LBGs in the HDF and concluded that 
these objects are dominated by young stellar populations 
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(<0.2 Gyr) and substantial amounts of dust (median 
E{B - y)~0.3). However, the bulk (11/17) of the ob- 
jects in their analysis came from the spectroscopic sam- 
ple of Lowenthal et al. (1997) who allowed objects that 
are redder — and so presumably more dusty — than 
those selected using the now so familiar criteria of Stei- 
del et al. that are used in our present KDF work. Indeed, 
Shapley et al. (2001) apphed the SED-fitting technique 
of Sawicki & Yee (1998) to a large sample of z~3 LBGs 
selected solely using the Steidel et al. selection criteria, 
and found a median E{B — V)—Q.16 (lower than Sawicki 
& Yee 1998) and concluded that LBGs undergo relatively 
short periods (50-100 Myr) of very intense star forma- 
tion followed by a more quiescent star-forming phase (see 
also Sawicki & Yee 1998). 
It is thus clear that at least bright {TZ<25) LBGs at 
are dominated by fairly short episodes of star forma- 
tion and significant dust obscuration. It remains unclear 
whether this is also the case at lower and higher red- 
shifts and at the fainter magnitudes that we reach in the 
KDF. Nevertheless, motivated by the results of Sawicki 
& Yee (1998) and Shapley et al. (2001), we take as our 
fiducial model the 100 Myr-old star-forming SED from 
the 1996 version of the Bruzual & Chariot (1993) and 
attenuate it with E{B - I^)=0.15 of dust. In § ITlUI 
we show that this fiducial model reproduces the Steidel 
et al. (1999, 2003, 2004) observed redshift distributions 
of z^2.2, z~3, and z~4 galaxy samples (the case for 
2'--^1.7 is less clear). Nevertheless, to monitor the impact 
of our choice of SED model on our LF results we carry 
out all our calculations in parallel, considering a grid of 
SED models that includes two stellar population ages, 
10 Myr and 100 Myr, and seven values of dust attenua- 
tion, E{B ~V)=0-0.3 in steps of 0.05. As we will discuss 
in 14.31 the dependence of the LF on these assumed dust 
and age values is negligibly small at z~4 and z~3, but 
becomes more significant at the lower redshifts. 

The sample completeness function p{m, z) 

Next, we must correct for incompleteness of our cat- 
alogs that is brought on both by the imperfect object 
detection efficiency and by the scatter of high-z galaxies 
across the boundaries of our color-color selection boxes. 
We do not (here or elsewhere in our LF calculation) 
explicitly correct for foreground, low-z interlopers that 
contaminate our high-2 sample. Such interloper contam- 
ination is known to be very small in the spectroscopic 
samples of Steidel et al. (1999, 2003, 2004). Because — 
as is discussed in KDF I — the contamination fraction 
is expected to at worst remain constant and in all likeli- 
hood to fall towards fainter magnitudes, the contamina- 
tion should be equally small or even smaller in our KDF 
data. 

We measure the amount of incompleteness by implant- 
ing simulated galaxies into our images and then seeking 
to recover them using the very same procedures that we 
used in making our data catalogs (see KDF I). Incom- 
pleteness (both detection incompleteness and the loss of 
galaxies due to scattering out of the high-z color-color se- 
lection boxes) is a function of apparent magnitude, with 
fainter galaxies suffering larger incompleteness than their 
brighter kin. It also depends on the true colors of our tar- 
get galaxies, and so their redshifts and intrinsic SEDs. 
We insert artificial objects with the expected colors of 



high-z galaxies generated as described in § 13.1.01 Our 
artificial objects have point-source profiles because — as 
we discussed in KDF I — the seeing in our images is suffi- 
ciently poor (FWHM'^l") to ensure that high-z galaxies 
are spatially unresolved. The recovered fractions form 
the completeness function p(m, z), which is the proba- 
bility that a galaxy of a given apparent magnitude (in / 
at z~ 4, TZ at z~3. Git at z~2.2, and G at z~1.7) and 
redshift z, dust attenuation, and model age, matches our 
sample selection criteria. 

The function p(m, z) is measured separately for each of 
the four redshift samples (z'--^4, z~3, z~2.2, and z^l.7) 
and — given the small differences in the image properties 
of our five KDF fields — is recalculated for each KDF 
field. The function is sampled in steps of Am=:0.5 in in- 
put apparent magnitude and Az=0.1 in redshift, and for 
the 8 combinations of age and reddening discussed above. 
At each step in this parameter grid several hundred sim- 
ulated objects are implanted at quasi-random locations 
in the image. These positions are always the same for the 
different steps in the parameter grid, but are otherwise 
unremarkable and sample the images fairly. 

We can use our calculated p{m, z) to test whether our 
assumptions about age and amount of dust are reason- 
able — i.e., whether we can reproduce the observed red- 
shift distributions of high-z populations. The shaded 
histograms in Fig. |21 show the redshift distributions of 
the spectroscopic samples of Steidel et al. (1999, 2003, 
2004). These Steidel et al. spectroscopic samples contain 
galaxies with a range of apparent magnitudes, and repre- 
sent the underlying color-selected population convolved 
with a spectroscopic success function that is not trivial 
to model. However, the bulk of their spectroscopic sam- 
ples consists of galaxies with 7?. ~25 — i.e., galaxies with 
photometric errors that are similar to those of TZ ^--^26. 5 
galaxies in the KDF. 

Thus, our p(m, z) for 7^=26.5 KDF galaxies should 
match the observed redshift distributions if our model- 
ing is a reasonable representation of reality. Fig |2] shows 
that this is indeed the case for our fiducial model with 
E{B - y)=0.15 and starburst age of 100 Myr (the 10 
Myr models, which are not shown, are also very good): 
Within each panel of Fig. [51 the solid curves show the 
function p(7?,=26.5, z) for the three values of E{B — V) 
(which generally increase from left to right), with the 
thicker curve marking E{B — T^)=0.15. Although our 
fiducial model does not give a uniquely-matching solu- 
tion (for example, a superposition of lower and higher 
E{B — V) values might work just as well), it does give a 
good agreement with the data. The fiducial model works 
remarkably well at z~4, 3, and 2.2; it works less well at 
Z'--^1.7 as it predicts a redshift distribution with a some- 
what lower median redshift than is observed, but so do 
all our other p(m, z) at that redshift (we will revisit this 
issue later). At any rate, the ability of our p(to, z) model- 
ing to reproduce the observed redshift distributions gives 
us confidence in that modeling, in the modeling of effec- 
tive survey volumes, Vef / , that are based on it in the next 
subsection, and thence in our estimate of the luminosity 
functions. 



The effective volumes 



Keck Deep Fields. II. The UV Luminosity Function 



5 



Finally, T4/ / is calculated for each of the five fields by 
integrating the probability function p(m, z) over redshift: 

dV 



dz 



■p{m, z)dz, 



(2) 



where dV/dz is the comoving volume per square ar- 
cminute in redshift slice dz at redshift z and ^/ is the 
area of the field in arcminutes. The Veff is calculated 
separately for each color-selected redshift sample (z~4, 
z~3, z~2.2, z^l.7), each galaxy model (which consist of 
the eight combinations of reddening and star burst age), 
each apparent magnitude step (Am=0.5 in J, TZ, GTZ, or 
G, depending on the redshift sample being considered), 
and each of the five KDF fields. 

3.2. The incompleteness- corrected number counts 

We next calculate the incompleteness-corrected galaxy 
density (f>(rn) as a function of apparent magnitude and for 
each Vef / model. We carry out the calculation separately 
for each of the redshift samples, z^A, z~3, 2:^2.2, and 
z^l.7. The calculation is first carried out independently 
for each KDF field and then the results are averaged 
together. 

For each field / we compute in 0.5-mag bins the num- 
ber of galaxies satisfying the color selection criteria, and 
then correct for incompleteness using the effective vol- 
ume: 

Here, Nf[m) is the number of observed galaxies within 
the magnitude bin m ± 0.25 in that field, T4//(m) is the 
effective volume of that field, cj) / (m) are the incomplete- 
ness corrected number counts in units of mag~^Mpc~'^ 
and the factor of 2 converts from counts in the 0.5 mag 
bins to counts per mag. 

The results of the individual fields are then weighted 
by field area, Af^ and averaged to yield the incomplete- 
ness-corrected galaxy number density for the entire KDF 
survey, 

- ^4^. (4) 

We restrict the calculation of the final KDF (t>{rn) to 
magnitude bins no fainter than the 50% detection com- 
pleteness limit in each field of the KDF. This limit is 
deeper by 0.5 mag for three of our fields than for the 
other two (see § |21 and KDF I) and so the average (/)(m) 
in Eq.^is computed using all five fields at all magnitudes 
except for the faintest magnitude bin in which only the 
deeper three fields are used. 

The uncertainty in <f>{m) combines two sources of un- 
certainty, namely the uncertainty in number statistics, 
S(f)]\[{m), and an estimate of field-to- field fluctuations, 
^4'f2f{iTi). The number statistics uncertainty is simply 
the Gaussian 54)M{jn) = (^(m)[7V(m)]~°-^, where N(m) is 
the total number of galaxies in all five (or three) fields 
in that magnitude bin. The field-to-field uncertainty 
5(f) f2 f is estimated using bootstrap resampling, whereby 
we generate 500 new realizations of (j>{m) via Eq. 0] but 
now in each realization choosing fields randomly with re- 
placement and so allowing the same field to be included 
more than once (or not at all) in a given realization. The 
field-to-field uncertainty (50/2/ is then taken to be the 



rms value of the 500 (j){m) resampled realizations. These 
two sources of uncertainty are then added in quadrature 
S(j) — [{S(j)N)'^ + (<5(/>/2/)^]°'^ to give us the total uncer- 
tainty. 

3.3. Absolute magnitudes 

The final step is to convert the (/)(to) into the luminos- 
ity function (j){M). As we discussed earlier, we compute 
the luminosity function at rest-frame 1700A both to min- 
imize fc-corrections, and to retain commonality with the 
work of Steidel et al. (1999) at brighter magnitudes. 

The absolute magnitude, M, is derived using the usual 
cosmological distance modulus, DM, and /c-correction, 
K, 

Afi7oo = - DM - K, (5) 

which we rewrite as 

Mi7oo = mA„,, - 5 log{DL/10pc) + 2.5 log(l + z) 

+ ("11700 - WA„i,,/(i+z)). (6) 

Here, Dl is the luminosity distance and rriAobs is the 
observed magnitude in the principal filter for the red- 
shift sample being considered (/, TZ, GTZ, and G for 
z~4, 3, 2.2, and 1.7, respectively). The last term of 
Eq. El (mi7oo — 'm\„tj(i+z))i is the /c-correction color 
between rest-frame 1700 A and the principal filter in the 
rest-frame for the redshift sample in question. This k- 
correction color is expected to be very small because of 
our decision to work at rest-frame 1700A. This expecta- 
tion is illustrated in the top two panels of Fig. O where we 
plot the /c-correction color for two representative galaxy 
models selected out of a larger ensemble that we tested. 
As Fig.Oillustrates, the value of the color term is indeed 
very close to zero for /-band at z~4, 7?.-band at z^3, the 
composite GTZ at z~2.2, and G-band at z~1.7. Over the 
redshift ranges selected by the color selection criteria we 
are using, the deviations from zero are typically no larger 
than 0.1 mag with a redshift-dependent Smag range of no 
more than ^±0.1 for a given model. These small offsets 
are negligible and so we set (toi7oo — ™Aob,/(i+-2)) ^"^^"^ 
in Eq. to arrive at 

Mi7oo = ™a„,, - 5 log(I?i/10pc) + 2.5 log(l -f z). (7) 

Applying Eq.Qto our (j){m) we at last arrive at the rest- 
frame 1700A luminosity function, (j){M). The resulting 
luminosity function measurements are described in the 
following section, § 0] 

4. THE OBSERVED LUMINOSITY FUNCTIONS 

4.1. Description of the luminosity functions 

The data points in Fig. ^ show our LF measured using 
our baseline Veff model of 100 Myr-old starburst with 
E{B — V^)=0.15 (dependence on model assumptions will 
be discussed in § 14.31) . In addition to the KDF data we 
also include the z~4 and z^3 LF points of Steidel et 
al. (1999). We stress that our analysis of the KDF data 
follows closely the selection and LF analysis procedures 
used by Steidel et al. (1999) for these brighter galaxies 
and so combining their results with ours should be ro- 
bust and free of systematic uncertainties. The KDF and 
the Steidel et al. samples complement each other: the 
Steidel et al. (1999) measurements provide good statis- 
tics at the bright end, but do not probe fainter than 
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-^^1700=— 20.5 at z~3 and MnQo——21 at z~4, while the 
KDF have good statistics at the faint end, probing ^1.5 
magnitudes deeper, but lack the statistics at the bright 
end. Unless otherwise stated, throughout the rest of this 
paper we will always combine the Steidel et al. (1999) 
measurements and our fainter KDF data when discussing 
the z~3 and z^-^A LPs. 

To give a more physical meaning to our luminosity 
scale we relate rest-frame UV magnitudes to star for- 
mation rates. Combining the relation between SFR and 
luminosity given by Kcnnicutt (1998; see also Madau, 
Pozzetti, & Dickinson 1998) with the definition of AB 
magnitudes (Oke 1974) gives 

SFR = 6.1 X 10-^^+°-^^'^'™°'>MQyr''\ (8) 

This conversion is valid in the absence of dust and for 
a stellar population that is forming stars continuously 
and whose UV light is dominated by massive, short-lived 
stars that are being produced on a Salpeter (1955) stel- 
lar initial mass function (IMP) with mass range 0.1< 
-^0^100. The SFR scales are plotted on the top axes in 
Pig. 0] The SFR scale should not be taken too literally 
because of the uncertainties in the assumptions underly- 
ing Eq.lHl Nevertheless, this scale gives us a useful, more 
physical reference frame for the luminosity functions. As 
Pig. El shows, the KDF LP reaches to Mi7oo=-18 at 
z~2.2 and Mi 700=— 19 at z~4. It thus rivals in depth 
LPs measured in the HDFs (e.g., Sawicki, Lin, Yee 1997; 
Steidel et al. 1999). Under the aforementioned assump- 
tions, our KDF LP reaches down to galaxies with star 
formation rates of ~2 Mq/jt and lower — comparable 
to the SFR in the Milky Way today. 

The solid curves in each panel of Pig.^jshow Schechter 
function fits to the data (Schechter 1976), and the gray 
shaded regions show the la uncertainties in those fits. 
We defer the description of the details of the Schechter 
function fitting to § 14.21 and here we use the fits only to 
guide the eye and to give a first comparison between LPs 
at different redshifts. 

Our LP is the best constrained of all the red- 
shift bins we consider and so, when comparing luminos- 
ity functions at different redshifts, we will use the z~3 
LP as reference. The dashed curves in the other panels 
of Fig^show the fit to this fiducial. Comparing the 
data at z^A with our fiducial z~3 curves immediately 
suggests that the LP undergoes evolution with redshift: 
The number density of faint galaxies at z~4 appears to 
be significantly lower than at z~3. At the same time, 
the number density of luminous galaxies app ears to re- 
main unchanged from z^A to z^3. In 14.31 and §|Slwe 
will explore in detail whether this LF evolution is real or 
simply the result of a selection effect or other artefacts 
(we conclude that it is very likely real at least between 
z~4 and z~3). 

For comparison. Pig. ^ also shows the recent low-z 
rest- 1500 A LP measurements from the GALEX mission. 
These LPs are shown as dotted curves in Fig. m with 
the left dotted curve showing the z~l GALEX LF and 
right one showing the z~0 LF (Arnouts et al. 2005 and 
Wyder et al. 2005, respectively). The strong evolution 
of the LP of star- forming galaxies from z~0 to z~l seen 
in the GALEX data has been recognized for some time 
(e.g., Lilly et al. 1995) and is responsible for the steep rise 
in the UV comoving luminosity density of the Universe 



over that redshift interval (e.g., Lilly et al. 1996; Schimi- 
novich et al. 2004). Similarly, the increase in the number 
of luminous galaxies from z~l to higher redshifts, z> 2 is 
part of the familiar UV-sketch of the cosmic star forma- 
tion history (e.g., Madau et al. 1996; Sawicki, Lin, & Yee 
1997; Giavalisco et al. 2004). In contrast to these well- 
known broad trends, the more subtle evolution in the LP 
from z^A through z~3 to z~2.2 that is revealed in our 
KDF data has not been so far explored because until now 
there was a lack of well-selected and well-tested samples 
that have good statistics over a wide range in luminosity. 

4.2. Parametric representation of the LF 

We fit the binned luminosity function data with the 
Schechter (1976) function, 

(bmodel{M)^^*${M) (9) 

= 0*O.4/n(lO)dex{[O.4(M* - A/)](i+")} 

which we evaluate over M*, 0*, and a. In practice, we 
do the fitting using minimization, where for a grid of 
M* and a values we compute a corresponding grid of 
values using 

x^(M^r,a) = g[ ^^-(^^yp'(^^Y ' (10) 

where the sum is taken over the M1700 magnitude bins, 
and (j)modei {M) are computed using Eq. Instead of 
adding a third, (j)* , dimension to the grid, the computa- 
tion is considerably accelerated by optimally calculating 
(f>* at each (M* , a) in the grid using the analytic relation 

EM^HM)/a^iM) ' 

which is derived by minimizing Eq. 1101 via dx^ /d(j)* = 
0. We then search the grid of values to select its 
minimum, Xmim adopt its corresponding M*, 0*, 
and a as the best-fit Schechter function parameters. The 
values of these best-fit parameters are listed in Tabled 
and plotted in Pig. [S] 

The error contours shown in Fig. |S1 are computed by 
recalculating the best-fitting M*, (fi* , and a, but now 
with sets of data values 4>data{M) that have been per- 
turbed randomly according to their standard deviations 
a{M). For each of the redshift bins, we generate 250 
such perturbed realizations and use their x^ to map out 
the regions of parameter space that correspond to the 
best-fitting 68.3% of such realizations. 

4.3. Systematic effects in the LF measurement 

Uncertainty due to k-corrections 

As we discussed in § 13.31 setting the fc-correction color 
term to zero introduces only a very small, ^0.1 mag, sys- 
tematic bias in the determination of absolute magnitude 
(see Pig. 13). Consequently, ^0.1 mag can be taken as 
the systematic uncertainty in our determination of the 
positions of 4>{M) bins. This uncertainty is too small to 
affect our LPs significantly and so we do not consider it 
further. 

The lack of spectroscopic redshift information for ob- 
jects in our sample may also introduce a systematic ef- 
fect, albeit — as we show here — a negligible one. Our 
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lack of spectroscopic redshifts means that we do not know 
whether a particular object in a given color-selected red- 
shift sample is near the lower or the higher end of the red- 
shift interval. This uncertainty may introduce a system- 
atic bias since at the same apparent magnitude intrinsi- 
cally less luminous objects are more likely to come from 
somewhat lower redshifts (and hence occupy a smaller 
effective volume) than intrinsically more luminous galax- 
ies. However, as is shown in the bottom panels of Fig.|21 
the effect of this redshift uncertainty on the value of the 
distance modulus DM and hence on the derived Afiyoo 
is small: it is ±0.2 mag for the and z^i samples, 
±0.4 mag at lower redshifts. Such systematic offsets of 
a few tenths of a magnitude in M1700 are comparable 
to uncertainties introduced into the LF measurement by 
N'^/'^ statistics and field-to-field variance as reflected, for 
example, in the uncertainties shown in Fig. |21 More- 
over, because the shapes of the LFs are similar in all 
redshift bins, the systematic bias introduced by this ef- 
fect works in the same direction at z~4, 3, 2.2, and 1.7, 
thereby reducing any systematic differences between the 
redshift bins. In summary, then, the biases introduced 
by our lack of spectroscopic redshifts are highly unlikely 
to drastically affect our results. 

Dependence on cosmology 

Throughout this paper we have assumed a cosmological 
model with (f^M, ^^A, i?o) = (0.3, 0.7, 70 km s'^ Mpc"!). 
However, it is well known that many quantities, including 
the luminosities and volume elements directly relevant 
in LF calculation, can depend strongly on the assumed 
cosmological model. Consequently, we feel it important 
to explore how the choice of cosmological model impacts 
our results. 

The choice of cosmology affects the derived luminos- 
ity function primarily through a change in the relative 
number density normalization (i.e., a change in (j); Eq.|^ 
and luminosity (i.e., a change in the absolute magnitude 
scale, M; Eq. [Tj). The dependence of the faint-end slope 
is weak as it enters only through a small difference in 
how the change in cosmology affects the derived Vef / at 
different apparent magnitudes. We have tested that the 
effect on the faint-end slope is negligible by recomputing 
our LFs for a model with (^m,^KtHq)={1, 0, 100 km 
s""'^ Mpc""'^) and consequently, we will focus our discus- 
sion on the effects on M and 

Figure shows how changing cosmology from our as- 
sumed (f7A/,f^A,ifo)=(0.3, 0.7, 70 km s^^ Mpc^i) to 
two other often-used cosmological models changes the 
resultant absolute magnitudes (top panel) and number 
densities (bottom panel). In both cases, the effect can 
be quite large: the difference in absolute magnitudes be- 
tween the two models can be as large as 1-2 mags, and 
the change in number density as large as an order of 
magnitude. However, the relative change in M and are 
only weakly dependent on redshift for a given cosmolog- 
ical model over the redshift range we study: relative to 
z~3, the z~1.7, z~2.2 and 2;~4 M scales are shifted by 
less than ^0.15 mag by a change of our assumed cos- 
mology; similarly, the number densities are also affected 
only weakly, at the <20% level. A direct consequence 
of this is that while the absolute (j) and M scales depend 
strongly on cosmology, the relative shapes and (j) and M 
normalizations from redshift bin to redshift bin are virtu- 



ally unchanged. As a result, any real evolutionary trends 
seen in the LF from to z^3 to z^2.2 to z^l.7 are 
virtually independent of the assumed cosmology. 

Dependence on galaxy model properties 

Our LF calculation relies — as it must — on our esti- 
mates of the effective volumes of the survey (see § OJ. 
These estimates are based on models of the ex- 

pected colors of high-z galaxies and therefore may de- 
pend critically on our model assumptions. Throughout 
this paper we have assumed a baseline model for high-z 
galaxies based on a 100-Myr old starburst reddened by 
E{B-V)=Q.lb of dust. As we have shown in S lTTlll this 
model reproduces well the observed redshift distributions 
of high-z galaxies, while models that are similar but have 
different amounts of dust do a poorer job. Nevertheless, 
we wish to explore how strongly model-dependent are 
our Veff estimates and the resultant LF measurements. 

To test how mo del- dependent our LF results are, we 
have repeated all the steps in our LF calculation (§ 13.11 — 
S I4.2|I for a grid of galaxy models that contains starbursts 
of two different ages (10 and 100 Myr) and eight values 
of E{B—V) (0-0.3 in steps of 0.05, and also a composite 
model that contains a mixture of galaxies with E{B — V) 
values drawn evenly from that E{B — V) range). We 
note that we have applied the different models to our 
own KDF points but were unable to do so to the points 
that come from Steidel et al. 1999 since we do not know 
in detail the Vef / model that was used in their LF work. 
Consequently, the effect of changing Vef / models will be 
seen only in our data and so will manifest itself fully at 
the faint end of the LFs, while the bright ends of the z~4 
and z~3 LFs, which are dominated by the Steidel et al. 
data, will remain unchanged. 

We show the effect of changing the Vef / model in two 
ways. First, the recomputed LFs are plotted in Fig.Q in 
which the data points at each M1700 show the results of 
7 LFs computed using the seven discrete E{B — V) val- 
ues, and the solid black curves show all the correspond- 
ing Schechter function fits; the gray curves in the z^4, 
2^^2.2, and z~1.7 panels reproduce the z~3 LF Schechter 
function fits for comparison. Second, the number density 
value of the Schechter function fits at a fixed magnitude, 
-^^1700=— 20, are shown as a function of E{B — V) in 
Fig.lHl here, the black lines connect the 100 Myr-old star- 
burst models and the lighter lines are for the 10 Myr-old 
ones. 

As can be seen in the z~3 panel of Fig. the choice 
of Veff galaxy model has a miniscule effect at z~3 — 
an effect of changing the Veff model has an effect that 
is no larger than the errors due to number statistics and 
cosmic variance (see Fig0J). Fig. |H1 confirms that at z^i 
the dependence on model is small: the number density of 
galaxies at Mi7oo=— 20 ranges over the range of models 
by <5% compared to our baseline model. We therefore 
conclude that at z~3 the LF we measure is very robust 
with respect to model assumptions. 

While the results at z~3 are virtually model- 
independent, the choice of Veff model has a larger (al- 
though still small) effect at z~4. Here, the 0m=-2o val- 
ues range by up to ±~20% with respect to our baseline 
model. However, despite this variation, the faint end of 
the z~4 LF always remains substantially below the z^3 
LF, and the variation between different models at z^A 
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is not larger than the uncertainties introduced by field- 
to-ficld variations and iV^/^ statistics, as can be seen by 
comparing Figs. Hand 01 We thus conclude that while 
the choice of Vef / model may possibly affect the details of 
the 2:~4 LF results, it is unlikely to alter the qualitative 
trends seen between the z^A and z~3 LFs. 

The choice of V^jf potentially has the largest effect 
at z~2.2. While models with E{B - V)>0 give a tight 
grouping of LF results, the two models with E{B — V)=0 
give LFs that have a substantially higher (/) normaliza- 
tion. These higher LFs are a direct result of the fact 
that the colors of z^2.2 models with E{B — V)=0 ap- 
proach very closely to the boundary of the color selec- 
tion box (see Figs. 4 and 5 in KDF I): large numbers of 
model galaxies in these models are scattered out of the 
color selection box, resulting in low p{m, z) values, low 
Veff, and hence high number densities (p in the resul- 
tant LFs. In our case t he eff ect is not critical because 
(as we have argued in S 13.1.0(1 the E{B — V)=G models 
are probably not realistic on other grounds; moreover, 
even the adoption of the extreme E{B — V)—Q results 
does negate the possibility that the evolution seen from 
z~4 to continues onwards to z~2.2 — in fact, an 
E{B — V^)=0 makes such an evolutionary trend stronger 
at z~2.2 than in the case of other E{B — V)\&\vies. How- 
ever, our E{B — V)=Q case at z^2.2 illustrates that in 
general one must be careful when computing luminosity 
functions of color-selected samples because such LFs can 
be very strongly dependent on the assumed properties of 
the high-z galaxy populations that one is trying to study. 

At zr^l.l the LF appears to have an even higher num- 
ber density normalization than at higher redshifts. We 
note however that we do not place much faith in our de- 
termination of the z^l.l LF as we discuss in more detail 
in S lOl 

In summary, we conclude that while the choice of Vef / 
models used in calculating the LF may have an effect on 
the LF results, in our case such effects are small: they are 
no larger than the uncertainties introduced by field-to- 
field variance and by \/7V statistics and in any case do not 
affect the qualitative evolutionary trends seen between 
z~4, 2^3. Meanwhile, the 2~2.2 model 14// may suffer 
from systematics that may adversely affect the accuracy 
of our LF measurements at those redshifts. 

Field-to-field fluctuations 

Galaxy clustering introduces field-to-field fluctuations 
in the galaxy distribution. This effect — often termed 
cosmic variance'^ — limits the accuracy with which the 
luminosity function can be measured. This is particularly 
true for surveys that rely on small single pointings, such 
as the HDF (e.g., Sawicki et al. 1997; Steidel et al. 1999), 
FDF (Gabasch et al. 2004) or the Hubble Ultra Deep 
Field (e.g., Bouwens et al. 2006). 

To mitigate the effects of cosmic variance, the KDF 
consists of five large fields grouped in three spatially- 
independent patches. When calculating the luminosity 

^ This name is not strictly correct as cosmic variance refers to 
the variance in samples that are fundamentally limited by the finite 
size of the Universe, such as, e.g., the largest-scale fluctuations 
on the Cosmic Microwave Background. Nevertheless, because it 
has entered the common nomenclature of galaxy evolution studies, 
we use the term cosmic variance interchangably with ficld-to-field 
variance throughout this paper. 



function we have taken field-to- field variance into account 
by including a bootstrap resampling measurement of this 
variance in our errorbars (see § 13. 2|) . Here, we illustrate 
the strength of field-to-field fluctuations by computing 
the characteristic galaxy number density 0* in each of 
the five KDF fields. We make the measurement of the 4>* 
while holding the shape of the luminosity function fixed. 
Specifically, at each redshift we hold a and M* constant 
at the values we measured earlier for the full dataset 
(Table ^1 while letting 0* be a free parameter. We also 
exclude the Steidel et al. (1999) z~4 and z~3 bright-end 
data from the measurement here. Table El summarizes 
the results and contrasts them with the (/)* values for the 
full dataset. The field-to- field fluctuations are generally 
small. They are largest at z~4 where the total number of 
objects in our sample is smallest: here, the largest excur- 
sion — that in the 03B field — is ^1.5 times the fiducial 
value, although the RMS scatter is significantly smaller 
than that and consistent with our bootstrap estimates. 

These field-to-field fluctuations, while relatively small 
in the KDF, underscore the need for multiple sightlines 
when determining the galaxy number density, luminosity 
function, and derived quantities such as luminosity and 
star formation rate densities of the universe. Single-field 
studies, especially if they rely on small fields such as the 
HDF or the UDF, have no way of monitoring this impor- 
tant source of error. In contrast, here, in the multi-field 
KDF, we have estimated field-to-field variance throught 
bootstrap resampling and included it explicitly in our 
error budget. 

4.4. Comparison with other surveys 

Over the last few years several authors have presented 
luminosity functions of UV-selected galaxies at high red- 
shift, z>\. The vast majority of these measurements 
were made using either full-blown photometric redshifts 
or the two-color selection techniques inspired by the suc- 
cess of the Steidel et al. surveys. The bulk of this work 
can be divided into two groups: those that use the very 
deep images of one or both Hubble Deep Fields (HDFs; 
Williams et al. 1996; Casertano et al. 2002) to probe the 
LF to very faint limits that are comparable to the KDF 
(e.g., Gwyn & Hartwick 1996; Sawicki et al. 1997; Stei- 
del et al. 1999, ), and those that use wider but shallower 
ground-based data (e.g., Steidel et al. 1999; Iwata et al. 
2003; Ouchi et al. 2004a; Gabasch et al. 2004). 

Here, we compare our 2;~3 and z~4 LFs with two re- 
cent wide but deep surveys, namely the z^A: LF of the 
Subaru Deep Survey (SDS; Ouchi et al. 2004a) and with 
the z~3 and z~4 LFs of the FORS Deep Field (FDF; 
Gabasch et al. 2004). We also compare our measurement 
of the LF's faint end with the Steidel et al. (1999) analysis 
of the HDF- North, for although the HDFs are too small 
to be truly adequate for LF determination, the Steidel et 
al. (1999) HDF-based faint-end slope of «=— 1.6 is often 
used in the literature, particularly by workers deriving 
the integrated UV luminosity density and star formation 
density of the Universe at higher redshifts (e.g., Bouwens 
et al. 2004; Bunker et al. 2004; Dickinson et al. 2004; Gi- 
avalisco et al. 2004; Yan & Windhorst 2004). 

The Hubble Deep Field LF of Steidel et al. (1999) 

Steidel et al. (1999) have applied a modified version of 
their color-color LBG selection criteria to the northern 
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Hubble Deep Field in order to extend to fainter magni- 
tudes their z~3 and z~4 LF analysis. Their z^3 and 
z^A HDF LFs are shown in Fig. as upward-pointing 
triangles. At z~3, our KDF LF (filled circles) is in good 
qualitative agreement with the Stcidcl ct al. 1999 HDF 
results albeit the KDF — with its much larger area and 
multiple pointings — provides a much more robust and 
precise measurement. Our faint-end slope Q!=— 1.43lQ'|!)g 
is significantly shallower than the a=— 1.6 reported by 
Steidel et al. from the HDF (or that by Sawicki et al. 1997 
who measured the high-z LFs in the HDF using a differ- 
ent selection technique). As we will discuss in Sawicki & 
Thompson (in prep., KDF III), the consequences of this 
difference in a are not negligible in calculations of the 
UV luminosity density and star formation density of the 
Universe. At z~4, the HDF data suffer from very poor 
statistics in addition to potential problems with cosmic 
variance, but they are lower than the z~3 points and, 
within their large errorbars, agree well with the KDF 
z-4 LF. 

The FORS Deep Field 

The FORS Deep Field reaches depths similar to those 
of the KDF in a single VLT pointing that covers 
-50 arcmin^, or 29% the area of the KDF. The fact that 
it comprises only a single field means that it is highly sub- 
ject to cosmic variance effects, but — unlike the HDFs 
— it is sufficiently large that it is not dominated by small 
number statistics at the faint end. Gabasch et al. (2004) 
have computed photometric redshifts in the FDF and 
used them to estimate the galaxy LF at several redshifts 
and rest-frame wavelenghts. They do not present a rest- 
frame 1700 A LF and so we compare our KDF results 
to their 1500 A LF. Both these rest- frame wavelenghts 
probe the part of galaxy SEDs dominated by young, mas- 
sive stars and so the systematic biases introduced by this 
slight difference in wavelenghts should not be large. 

Figure pluses downward-pointing red triangles to show 
the 1500A Gabash et al. (2004) FDF luminosity functions 
at z~3 and z~4. The agreement between the KDF and 
FDF is very good, particularly when one considers that 
the FDF points are based on a single moderately-sized 
field and so do not account for cosmic variance and, more- 
over, that the FDF and KDF use very different means 
of sample selection. In particular, we note that both 
the FDF and the KDF show relatively shallow faint-end 
slopes at both z— 3 and z— 4 although we feel that the 
KDF result is much more robust both with respect to 
cosmic variance (because of its larger area and multiple 
independent fields) and sample selection (because of the 
extensive spectroscopic verification of its selection tech- 
nique). 

The Subaru Deep Survey 

The Subaru Deep Survey (SDS) consists of two large 
fields (the Subaru Deep Field, SDF, and the Subaru 
XMM Deep Field, SXDF). Both fields have been imaged 
with the Suprimecam large mosaic imager (600 arcmin^ 
per field) on Subaru and contain statistically large num- 
bers of galaxies and are very unlikely to be affected by 
cosmic variance. The SDS BVRi' images allow the se- 
lection of z— 4 and z~5 LBG samples and Ouchi et al. 
(2004a) have used such selection to estimate the LBG 



LFs at these redshifts. A strong limitation of the SDS 
selection of high-z galaxies is that their color-color selec- 
tion technique is not calibrated spectroscopically (only a 
small handful of z>l spectroscopic redshifts is known in 
the SDS fields) and must instead rely on models of galaxy 
colors for its definition. In Fig. [HI we use blue squares to 
plot the z~4 rest-frame 1700 A luminosity function from 
the SDF field (we omit the SXDF field which shows vir- 
tually identical results but is ~0.5 mag shallower than 
the SDF). At the bright end, the SDS z-4 LF agrees 
both with the KDF luminisity function (which at these 
magnitudes is dominated by the Steidel et al. 1999 data) 
and with the FDF result. However, at the faint end, the 
SDS has a much steeper LF slope than either the KDF 
or the FDF. The origin of this discrepancy is not imme- 
diately clear. One possible explanation is that at faint 
magnitudes, where as we have seen in S l4.3.0l the scatter 
in galaxy photometry can strongly affect the estimated 
effective volumes used in calculating the LF. Ouchi et al. 
(2004a) do not discuss what galaxy models they used to 
estimate their 14//, but it is possible that while our T4// 
estimates at z~4 are robust to changes in the assumed 
galaxy SEDs, their effective volumes computed for their 
BRi' color-color selection are less so. If this is the case, 
then the Ve// in the SDS work may well be underesti- 
mated resulting in overcorrections to the SDS LF at the 
fainter magnitudes. 

5. EVOLUTION OF THE LUMINOSITY FUNCTION 

In this section we compare our LFs at different red- 
shifts to search for signs of evolution. In what follows, 
we concentrate on Figs. EH and which compare di- 
rectly our z~4, 3, 2.2, 1.7 LFs that were determined us- 
ing the baseline 14// model (lOOMyr-old starburst with 
E{B — V)=0.15). As we will elaborate below, our z~3 
and z^4 LFs are the most unambiguous and free of sys- 
tematics, whereas the z~2.2 and especially z~1.7 LFs 
are likely biased; for this reason in Fig. 1101 we present 
the z~4 and z~3 LFs only, without the z~2.2 and z~1.7 
LFs that are shown in Fig. EH Thus, Fig. 1 101 shows only 
the most solid, bias- free results, while Fig.^Jshould thus 
be regarded as presenting all the results, including those 
biased by systematic effects. 

Our LF is most unambiguously constrained at 2^3, 
where the combination of KDF and Steidel et al. (1999) 
data covers the largest range in luminosity, is least de- 
pendant on the details of 14// modeling (see § I4.3.0|l . 
and puts the tightest bounds on the Schechter function 
parameters. Consequently, we will use the z~3 LF as our 
fiducial reference and compare the other redshift bins to 
it. 

The top panels of Figs.llllandllOloverplot the data and 
the Schechter function fits, and the bottom panels fur- 
ther highlights the differences between the three redshift 
samples by showing the data after they have been di- 
vided by the Schechter function fit to the z~3 LF. Three 
interesting evolutionary effects can be seen in Figs. ITTI 
andlTUl 

• A strong increase in the number density of low- 
luminosity LBGs from z^4 to z^3. 

• An accompanying apparent lack of change in the 
number density of luminous LBGs from z~4 to 
z~3. 
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• A possible continuation of the increase in the num- 
ber density of low luminosity galaxies from z~3 to 
z~2.2 and z-1.7. 

We discuss these three effects in the following subsec- 
tions. In each subsection, we give a brief phcnomcno- 
logical description of the effect and then concentrate on 
examining whether the evolution may be explained as an 
observational effect or whether it instead reflects a true 
evolution of the underlying galaxy population. Discus- 
sion of the possible implications for galaxy evolution of 
the observed LF evolution are deferred to § 

5.1. Increase in the number of low-luminosity galaxies 
from z^4 to z^3 

The most striking effect in Fig. ^| is the increase in 
the number density of low-luminosity galaxies from z^4 
to z~3. For our adopted fiducial E{B — V)=Q.lh model, 
the number of z^4: galaxies with — 21<Mi7oo^— 18 (i.e, 
within M*<Mi7oo<Af*+2) is only 0.44±0.05 of the 
number of such galaxies at z~3. In other words, the 
number of low luminosity galaxies increases by a factor 
of 2.3 from z~4 to z~3 — a result that is statistically 
significant at the ^llcr level. 

Is this increase a reflection of true galaxy evolution, or 
is it simply an artefact of some observational bias? 

Can cosmic variance be responsible? 

The effects of large scale structure can be a problem in 
small-area surveys such as the HDFs or those that con- 
sist of a single pointing, such as the FORS Deep Field. 
The KDF, however, consists of three widely-separated 
patches on the sky that probe statistically independent 
parts of the Universe. As was discussed in §§13 and 
01 field-to-field fluctuations are relatively small and the 
errorbars in our LF measurement include a bootstrap 
estimate of field-to-field variance of our survey. Given 
that the increase in the number of low-luminosity galax- 
ies from zr^A to z~3 is far outside these errorbars, we 
conclude that cosmic variance is unlikely to be respon- 
sible for the observed evolution of the faint end of the 
LF. 

Is our modeling of Vef / responsible ? 

As we discussed in § 14.3.01 the computation of the LF 
is dependent on the details of the Vef f mod eling of the 
survey. However, as we have shown in S 14. 3. 01 at zr^A and 
z~3 the change in the LF is remarkably small under a 
wide range of reasonable assumptions. As is illustrated 
in the top panel of Fig. [7\ varying the Vef / model as- 
sumptions cannot bring the faint ends of the z^3 and 
z~4 LF into agreement. And as Fig. |H1 shows further, 
it appears to be impossible to bring the number density 
of ^M* + l galaxies at z^3 and z^A into agreement by 
varying Veff model parameters: even in the most ex- 
treme E{B — V)=0 case, where the number densities at 
z~4 and 3 are closest to each other, there remains a very 
significant deficit of faint galaxies at z~4. We therefore 
conclude that the evolution of the faint end seen in our 
data is unlikely to be an artefact of the assumptions that 
underlie our calculation of Vef / . 

Can differential sample selection be responsible? 



The color selection criteria used to select galaxies at 
z~3 and z~4 are designed to select galaxies with simi- 
lar underlying SEDs. However, the color-color selection 
regions are not identical in relation to intrinsic galaxy 
colors (see Fig. and so we may question whether the 
deficit of sub-i* galaxies at z^4 results simply because 
at z~4 we are missing galaxies whose SEDs are such that 
they would have been included in the z~3 selection cri- 
teria. 

If there are significant numbers of such galaxies missing 
from our sample at z^4 then we should be able to include 
them by expanding the z~4 color-color selection region. 
However, the deficit at the faint end is so large — a 
factor of 2.3 — that no reasonable adjustment to the 
z^A color selection criteria can remedy it. As Fig. ^ 
shows, the galaxy color model tracks allow the possibility 
of modifying the z~4 color-color selection region by a few 
tenths of a magnitude. However, we have checked that 
an increase of even as much as 0.5 mag in both G~TZ and 
TZ — I would not be sufficient to bring enough galaxies 
into the sample to make up the factor-of-2.3 deficit (see 
also Fig. 4 in KDF I). Moreover, even if such a large 
change to the color selection criteria were permissible, 
it would result in an automatic increase of the effective 
volume Vef / that would largely counteract any gain from 
the increase in galaxy numbers. 

A further, secondary argument against differential 
sample selection relies on the fact that the evolution from 
z~4 to Z'^S appears to be differential with luminosity 
(see § I5.2|l . Since at a given redshift bright and faint 
LBGs are selected in the same way, a sample selection 
bias should result in a similar deficit of bright galaxies 
as of faint ones at z~4. The fact that no strong deficit 
is seen at the bright end of the z~4 LF further strength- 
ens the case that the evolution of the faint end is not an 
artefact of sample selection but is due to real differences 
between luminous and low-luminosity galaxies. 

Overall, we believe that the deficit of faint galaxies at 
z~4 is too large to be accounted for by differences in sam- 
ple selection and is most likely a real effect that reflects 
an underlying evolutionary change in the population of 
sub-i* galaxies over the ^600 Myr from z~4 to z~3. We 
conclude that the observed evolution of the faint end of 
the LF from z~4 to z~3 is likely real. 

5.2. Luminosity-dependent evolution from Z'^^ to z^3 

Figure El suggests that while there is a deficit of z~4 
sub-L* galaxies, the number of luminous galaxies re- 
mains virtually unchanged from Z'^4 to z~3. If real, 
this differential, luminosity-dependent evolution of the 
LF hints at important differences in how galaxies of dif- 
ferent luminosity evolve at high redshift. We examine 
the potential implications of this differential evolution in 
§0 but first, here, we ask if this differential effect is in 
fact real. 

Is the effect statistically significant? 

While there can be little doubt of the deficit of faint 
galaxies at z~4, the situation at the bright end is less 
clear because of the larger uncertainties in the individ- 
ual data points. In view of this problem, it is tempting 
to rely on the Schechter functions to compare the bright- 
end LFs, and indeed the Z'^^ and 2~4 Schechter function 
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fits are in excellent agreement at the bright end. How- 
ever, such a comparison can be misleading because the 
Schechter function fits incorporate data at all luminosi- 
ties and so at the bright end the fits may well be biased 
by the weight of the faint-end data where the statistics 
are so much better. It would be more robust to compare 
the z~3 and Z'^A number densities directly, but this is 
not straightforward because the data at z~4 and z~3 are 
sampled differently and cover somewhat different magni- 
tude ranges in the two redshift bins. 

To overcome these difficulties, we instead compare how 
the data at z~3 and 2:~4 deviate from our Schechter 
function fit. We proceed as follows. First, we calculate 
the quantity <I>(M), which is the ratio between the data 
and the 2^3 Schechter function fit, 

$(M) = (l)aata{M)/cP)ZHM). (12) 

The bottom panels of Figs . ITTI and [TUl show <i>(M) for our 
four redshift samples. Note that $(Af) is computed for 
each of the four redshift samples, but in each case the 
data are divided by the same, z^3, Schechter function. 
Consequently, $(Af) is close to 1 at z~3 (reflecting the 
fact that the 2;~3 Schechter function is a good represen- 
tation of the 2~3 data) but deviates from 1 for the other 
redshift bins, and in particular for the sample. We 
then compute the average $(Af), namely $, for galaxies 
brighter and fainter than Mi7oo=— 21.0 (i.e., ~M*) at 
both redshifts of_interest here, and z^A. Finally, 
the ratio of the $ at the two redshifts then tells us the 
amount of number density evolution that the given pop- 
ulation undergoes. By comparing the ratios of the $, we 
effectively cancel out the dependence of our comparison 
on the Z'^S Schechter function fit and are comparing the 
data at Z'^S and z^A directly. 

Figure IT^ shows these $ ratios for the bright (horizon- 
tal axis) and faint (vertical axis) ends of the luminosity 
function. The quantities in Fig. 1121 are always shown 
as evolution with respect to the z~3 case (i.e., we plot 
^z/^zr^s)- Three evolutionary scenarios are marked for 
reference: locations on the vertical straight line indicate 
no number density change in the bright end of the LF, lo- 
cations on the horizontal solid line indicate no evolution 
in the faint end, and locations on the diagonal line in- 
dicate equal number density evolution at the bright and 
faint end. The intersection of the three lines at (1, 1) 
marks the case of a non-evolving LF. 

As Fig. ^] shows, there is substantial change in the 
iiumber density ratio of sub-L* galaxies from z~4 to z~3: 
$(z~4) /$(z--3)=0.44±0.05 for faint galaxies, indicating 
that there is a 2.3-fold increase in the number density of 
faint galaxies that is statistically significant at the ^lla 
level (this is the faint-end evolution we discussed in §|OJ. 
At the same time, however. Fig. E| also shows that the 
number density of luminous galaxies at z^4 is virtually 
unchanged with $(z~4)/$(z~3)=0.87±0.19. 

To properly test whether the evolution of the LF is 
differential with luminosity, we must of course consider 
the joint uncertainty for the bright- and faint-end cases. 
The joint la and 2a uncertainties are illustrated by the 
error ellipses in Fig.^l and in this context, the distance 
of the z~4— >3 point from the diagonal "equal evolution" 
line indicates the amount of differential evolution. We 
find that 98.5% random realizations of the data in Fig. 1121 
are above the diagonal "equal evolution" line and so in- 



consistent with the differential evolution scenario. This 
strongly suggests that the high-z galaxy population is 
undergoing differential, luminosity-dependent evolution. 
Stronger confirmation of this assertion will require im- 
proved constraints on the bright end of the LF; this con- 
firmation will require LBG surveys with areas of several 
square degrees — i.e. an order of magnitude greater than 
used in the work of Steidel et al. (1999) that provides the 
bulk of the statistics at the bright end of the LF here. 

Can differential sample selection be responsible? 

The bright and faint galaxies in a given redshift bin 
are selected using identical selection criteria. Moreover, 
given that the bright-end galaxies come primarily from 
the Steidel et al. (1999) work and the faint-end ones ex- 
clusively from the significantly deeper KDF, they suffer 
from similar amounts of photometric scattering out of 
the color-color selection regions (at any rate, this scatter 
is accounted for through the Ve// approach). Conse- 
quently, as long as the bright and faint samples within 
that bin consist of the same mix of galaxy SEDs, any 
systematic effect should be reflected in the same way at 
both the bright and faint ends. 

At the bright end, the LF is constrained mainly by 
the Steidel et al. (1999) results, while the faint end is 
dominated by the KDF data. It is possible that despite 
our great care some subtle, unknown selection effect re- 
mains between the two samples. However, we consider 
this possibility to be extremely unlikely in view of the 
virtually identical observational approaches and sample 
selection and analysis techniques that were applied to the 
two datasets. 

Finally, we cannot rule out the possibility that the faint 
and bright galaxies at a given redshift represent pop- 
ulations that differ not just in luminosity but in some 
other property such as, for example, the amount of inter- 
stellar dust. However, if such true, intrinsic luminosity- 
dependent differences do exist, they do not alter our con- 
clusion that the LF evolution is differential with lumi- 
nosity, but merely shift the focus to a more specific, but 
still real, evolutionary effect. We conclude that while 
the different composition of the samples with luminos- 
ity within a given redshift bin may explain the differen- 
tial, luminosity-dependent evolution we see from z^4 to 
z~3, such differences in intrinsic galaxy properties would 
only underscore the point that the luminous and faint 
UV-selected galaxy sub-populations are different and not 
simply scaled analogs of each other. 

5.3. Evolution from z^3 to z^2.2 

At z~2.2 (and z^l.7) we do not have sufficient statis- 
tics in the KDF data to say much about the LF for galax- 
ies brighter than L*, nor — unlike at z~3 and z~4 — 
are there published results from the Steidel group that we 
could use to extend our luminosity range. Consequently, 
we focus on the evolution of the faint end of the LF only. 

FigurelTTlshows evolution of the LF from z^3 to z^2.2 
(see also Figs. 0] and Ej). The LF appears to undergo only 
a small amount of evolution from z~3 to z~2.2 for our 
fiducial Veff model (100 Myr-old starburst with E{B — 
F)=0.15): there is only a factor of l.liO.l increase in the 
number density of sub-L* with cosmic time. However, as 
Fig.lHlillustrates, our chosen fiducial Veff model happens 
to give nearly the minimal evolution in the LF from Z'^3 
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to z~2.2 and applying z~2.2 Ve// models with either 
larger or smaller extinction values results in stronger LF 
evolution. For example, using the E{B — V)=0.3 model 
for both z~3 and 2;~2.2 gives a factor of ^1.2 increase in 
the number density of sub-L* galaxies and a yet stronger 
evolution (factor of 2.3 number density increase) results 
with the E{B — V)=0 model. Likewise, adopting a Vgff 
model which assumes a galaxy population with a range 
of E{B — V) values gives evolution that is stronger than 
our fiducial case: using a flat number distribution over 
E{B~V)^0-0.3 results in an increase of 1.24±0.08 in the 
number density of galaxies fainter than Mi7oo=— 21.0. 

The above examples illustrate that the measurement 
of the z~2.2 luminosity function has a systematic un- 
certainty that depends on our assumptions about the 
properties of the UV-bright galaxy population at these 
redshifts. Whereas the and, especially, z~3 LFs are 
only weakly dependent on the assumed Veff model, the 
z~2.2 LF shows a much stronger dependence. A direct 
result of this dependence is that we cannot unambigu- 
ously determine the z^2.2 LF and the amount of LF 
evolution from z~3 to z~2.2. 

Despite these limitations, we can still put some useful 
constraints on the evolution of the LF over the 800 Myr 
from z~3 to z~2.2. As Fig. [3 shows (see also Fig. ISJ, 
there are at least as many galaxies at a given UV lumi- 
nosity at z~2.2 as there are at z~3. Depending on the 
adopted Ve// models, the number density can be close 
to identical at the two redshifts, but potentially can be 
up to a factor of ^2 higher at z^2.2 than at z~3. Note 
that this conclusion could potentially be further affected 
by differential sample selection between z~3 and z~2.2, 
although it is unlikely that sample selection differences 
will drastically modify the results given that the z~3 and 
z~2.2 color-color selection criteria select galaxies with 
similar intrinsic SEDs (Fig. ^ and are likely to capture 
the bulk of UV-bright galaxies at these redshifts, as can 
be seen in Fig. 5 of KDF I. 

In summary, we can conclude that there are probably 
at least as many — and possibly more — sub-L* UV- 
bright galaxies at z~2.2 as there are at z~3. There is 
no evidence for a decline in galaxy numbers with cosmic 
time. 

5.4. Towards 1.7? 

Figure shows apparent strong evolution of the LF 
from higher redshifts down to z~1.7. However, we feel 
that the LF measurement at z^l.7 is not reliable given 
systematic uncertainties in our LF estimation at this red- 
shift. The main systematic problem stems from the lack 
of robustness in the modeling of V^ff for the z~1.7 sam- 
ple. While we have high confidence that our V^ff mod- 
eling is robust at Z'--^4 and z~3, and somewhat less accu- 
rate but still partially reliable at z^2.2, we have far less 
confidence in its reliability at Z'-^^l.?, as explained below. 

While our Ve// modeling for the higher redshift bins 
accurately reproduces the observed redshift distribution 
of galaxies at z~4, 3, and 2.2, it fails to do so at z~1.7, 
as can be seen in Fig. |21 Another manifestation of this 
problem can be seen in Fig. ^ where the model colors 
of z=1.7 galaxies (z=1.7 is the mean redshift of the ob- 
served spectroscopic sample) are too red in [/„ — G com- 
pared to the color-color selection region defined by Stei- 
del et al. 2004. The upshot of these mismatches is that we 



do not have confidence that our modeling of V/ / accu- 
rately reflects the true volume from which the color-color 
selected galaxies in the z~1.7 sample are drawn. An inac- 
curate calculation of the survey volume has the potential 
to strongly affect the computed luminosity function. If, 
for example, our Ve// modeling underestimates the vol- 
ume of the survey, then this will translate into an over- 
estimate of the number density of galaxies. This may 
well be happening and would explain the high density 
normalization of the z~1.7 LF seen in Fig. 1111 

There are several possible underlying reasons why our 
Veff modeling at z~1.7 may be inaccurate. We do not 
at present know what E{B — V) value or star forma- 
tion history are appropriate for z~1.7 galaxies. At z^3, 
good estimates of these quantities are known from obser- 
vations (see S I3.1.0|) . and it's not unreasonable that sim- 
ilar E{B — V) and age values are also applicable at z^4 
and Z'--^2.2 given that these redshift bins are only 0.6 Gyr 
and 0.8 Gyr away from z~3. By z~1.7, however, 1.8 Gyr 
after z~3, the galaxy population may have significantly 
different properties than it does at z~3 as galaxies evolve 
towards the more prosaic, less star-burst dominated and 
less dust-obscured galaxy population observed by Z'^l. 
Assuming the Z'^Z reddening and starburst age values for 
the z^l.7 population may well strongly skew our Ve/ / es- 
timate and be refiected in the mismatches seen at z~1.7 
in Figs. ID and El 

In addition to systematic problems with the Ve// cal- 
culation, a second potential systematic effect may be 
affecting our LF measurement at z~1.7: it is possible 
that our z~1.7 color-color selected sample may be sig- 
nificantly contaminated by low-z galaxies. While the 
spectroscopically-determined contamination rate in the 
7^~24-25.5 z-1.7 sample of Steidel et al. (2004) is less 
than 5%, even small systematic offsets in the color-scale 
calibration between our and their photometric systems 
can result in a drastically increased contamination frac- 
tion. We stress here that our modeling of Ve// accounts 
for the scatter of high-z galaxies between the different 
high-z color-color selection regions, but does not account 
for the scatter of low-redshift (z<l) galaxies into the 
high-z selection boxes. At z^4, 3, and 2.2, the scatter 
from low redshift is unlikely to be an issue given that 
these selection regions are far- removed from the z<l loci 
in color-color space. In contrast, the Z'^1.7 selection re- 
gion is in close proximity to the region of space domi- 
nated by z<l galaxies and it is possible that the low-z 
interloper contamination rate may be quite high. A high 
foreground contamination fraction would result directly 
in an overestimate of the number density of z~1.7 galax- 
ies and an inaccurate luminosity function. 

In summary, unlike at Z'--^4, 3, and 2.2., we are not 
confident in our determination of the luminosity function 
at Z'--^1.7. We have chosen to present the Z'--^1.7 LF here 
for completeness, but we urge the reader to regard it with 
much caution. 

6. DISCUSSION: THE EVOLVING GALAXY POPULATION 

6.1. Summary of the observational evidence 

To summarize §[3 the number of faint (sub-L*) galax- 
ies increases with cosmic time, with ^2.3 times more 
galaxies of a given luminosity at z^3 than at z^4; this 
increase is statistically significant at the llcr level and 
— as we have argued in § 15.11 — it seems very unlikely 



Keck Deep Fields. II. The UV Luminosity Function 



13 



that it is an artefact of differential sample selection or 
cosmic variance in our data. At the same time, the evo- 
lution from 2~4 to z~3 appears to be differential with 
luminosity, since while the faint end of the LF evolves 
significantly the bright end appears to remain virtually 
unchanged. This differential effect is statistically signifi- 
cant at the 97% level and is unlikely to be due to a sys- 
tematic bias. The case for evolution from to lower 
redshifts is less clear because of potential systematic bi- 
ases, although it is possible that the number density of 
sub-i* galaxies continues to increase to at least 2~2.2; 
our analysis suggests that it is extremely unlikely that 
the number density of sub-L* galaxies at z~2.2 is lower 
than at 

Because our results at and Z'^A are highly robust 
(while those at lower redshifts are less so), in what follows 
we focus on the evolution of these two epochs. 

6.2. The importance of luminosity function evolution 

The evolution of the luminosity function measures only 
the evolution of the galaxy population as a whole and 
does not necessarily imply a direct correspondence in the 
evolution of individual galaxies. The observed evolution 
of the LP's faint end can be interpreted equally well as a 
change in the number density or luminosity of the observ- 
able population. However, individual galaxies are free to 
change their luminosities following trajectories that are 
far more complicated than a direct increase in luminosity 
or number density as the evolution of the faint end of the 
LF might naively suggest. Likewise, the apparent con- 
stancy of the LF's bright end does not necessarily imply 
that the luminous galaxies that populate it do not them- 
selves evolve. Evolution of the LF is clearly not a direct 
probe of the evolution of individual galaxies. 

However, while the luminosity or number density evo- 
lution of the LF does not necessarily reflect a direct 
corresponding evolution in the properties of individual 
galaxies, the fact that the LF does evolve constitutes 
an important suggestion that its constituent individual 
galaxies do evolve over time. Furthermore, the fact that 
the LF's evolution appears to be luminosity-dependent 
suggests that the evolution of individual galaxies is also 
differential with luminosity. It suggests that there may 
be real physical differences between low-luminosity and 
high-luminosity systems in properties such as the supply 
of gas available for star formation, the merger rates that 
may trigger such star formation, properties of the dust 
that obscures it, or the effectiveness of feedback that can 
regulate it. 

6.3. Some evolutionary speculations 

Both semi-analytic and SPH galaxy evolution models 
have been used to predict the shape of the LBG luminos- 
ity function at different redshifts (e.g., Somerville, Pri- 
mack, & Faber 2001; Nagamine et al. 2004). However, 
such models tend to produce a relatively constant, un- 
evolving luminosity function — we are not aware of any 
predictions in the literature for the evolution of the sub- 
L* end of the LF that we observe from to z^3. It 
would be interesting to see what modifications to these 
sophisticated galaxy formation models can reproduce a 
differentially evolving LF. In the meantime, in the ab- 
sence of such predictions, we turn to some simple phe- 
nomenological speculations about the possible nature of 



the evolution of individual galaxies that underlies the 
observed evolution in the LF. We explore three heuris- 
tic evolutionary scenarios which we use to illustrate how 
varying some simple properties of individual galexies can 
mirror the observed differential, luminosity-dependent 
evolution of the LF. 

Our three heuristic models are motivated as follows. 
SED studies suggest that high-z galaxies likely undergo 
episodes of intense star formation followed by more qui- 
escent periods and it is plausible that such episodes oc- 
cur several times in the life of a high-z galaxy (sec, e.g., 
Sawicki, & Yee 1998; Shapley et al. 2001). Moreover, 
it is clear that UV-sclccted liigh-z galaxies are obscin-ed 
by significant amoimts of starlight-absorbing dust (e.g., 
Meurer et al. 1997; Sawicki & Yee 1998; Ouchi et al. 
1999; Shapley et al. 2001; Papovich, Dickinson, & Fer- 
guson 2001; Vijh, Witt, & Gordon 2003) and there is no 
reason to believe that the properties of this dust such 
as its opacity or large-scale geometry - remain constant 
with time or star formation rate. Two of our heuris- 
tic scenarios (Scenarios B and C) are thus motivated by 
the possibility that the properties of dust or of the star- 
bursting episodes evolve with time and/or luminosity. 
Meanwhile, Scenario A investigates the more simple pic- 
ture that in the evolution of the faint end we are seeing 
the very first appearance of many of the low-luminosity 
galaxies. 

We stress that our three heuristic scenarios are not 
meant as a comprehensive survey of all possible evolu- 
tionary mechanisms. Clearly there are many others, but 
we focus on these three to illustrate some interesting pos- 
sibilities and motivate future foUowup studies. 

Scenario A: The first appearance of low-luminosity 

galaxies ? 

One of the simplest possible pictures of Lyman Break 
Galaxies is that they are objects that form their stars at 
a constant, unvarying rate for long periods of time. SED 
modeling can be used to constrain the ages of the ongoing 
episodes of star formation in LBGs and while such ages 
are notoriously dependent on the assumed star formation 
histories, it is the constant star formation rate scenarios 
that yield the oldest ages (Sawicki & Yee 1998; Papovich 
et al. 2001) thus providing upper bounds. Shapley et al. 
(2001) have modeled a large sample of relatively lumi- 
nous (typically 7?.~24) 2;~3 LBGs under the assumption 
of constant star formation. Of 72 luminous z^3 galaxies 
in their analysis, only 38% have ages older than 0.6 Gyr 
and so only 38% of the luminous LBGs seen at z^3 were 
forming stars at z^A — the remainder must have begun 
their current episodes of star formation more recently 
than z^4. Meanwhile, since the number density of lumi- 
nous LBGs is the same at z~4 as at z^3, 62% of lumi- 
nous 2~4 LBGs must have ceased star formation before 
Z'^3 to be replaced by the younger starbursts. Similar 
reasoning applied by Iwata et al. (2003) to an even ear- 
lier epoch suggests that only 20% of luminous LBGs at 
z^5 can still be seen at z^3 while the remaining 80% 
must have been replaced to keep the bright end of the 
luminosity function constant. These arguments suggest 
that individual LBGs cannot be in a steady star-forming 
state but at best are in a quasi-steady state where indi- 
vidual galaxies fade in and out of a given magnitude bin 
to keep the number density constant, even if timescales 
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for such fading are fong. Indeed, a detailed analysis by 
Ferguson, Dickinson, & Papovich (2002) calculated that 
the star formation histories of (luminous) z~3 LBGs are 
inconsistent with the observed number density of these 
galaxies at 2^4 unless episodic bursts of star formation 
are invoked. 

However, the situation may well be different for low 
luminosity galaxies — there it is possible to reproduce 
the observed LF evolution by postulating that a low- 
luminosity LEG remains at a constant UV luminosity 
once it start forming stars. The age distribution of sub- 
L* LBGs has not yet been measured, but let us as- 
sume that — as for the luminous LBGs — only 38% of 
low-luminosity z~3 LBGs are old enough to have been 
present at 2;~4. If at the same time we assume that 
none of the z^A LBGs have ceased star formation to fade 
out of the sample, then we can reproduce the ~2. 3-fold 
number density evolution of low-luminosity LBGs by as- 
suming that new low-luminosity LBGs are being simply 
added to the population between 2;~4 and z~3. 

Under this scenario it is possible that we are seeing 
large numbers of low-luminosity LBGs "light-up" for the 
first time in the epoch between 2;~4 and z~3. If this is 
the case then we may expect them to have low metal- 
licities — a property that should allow us to test this 
scenario through future observations. This picture is 
perhaps somewhat akin to "donwsizing" scenarios which 
postulate that star-formation activity shifts to lower- 
luminosity objects over time. It is possible, however, to 
explain the observed LF evolution with other evolution- 
ary models as we discuss next. 

Scenario B: Evolution in the properties of starforming 
bursts ? 

The LF evolution may also (or instead) be related 
to the frequecies, durations, or intensities of the star- 
bursting episodes that likely rule high-z galaxies. While 
the constant star formation assumption in SED fitting 
of LBGs can result in relatively old starburst ages, other 
assumed star formation histories can yield significantly 
shorter star-formation episodes (Sawicki & Yee 1998; Pa- 
povich et al 2001) and it is plausible that such periods of 
relatively brief elevated star formation may reccur several 
times in each galaxy between z~4 and Such fluctu- 
ating star formation rates are a feature of some models 
of Lyman Break Galaxies (e.g., Nagamine et al. 2004). 
As the SFR in a galaxy fluctuates over time, that galaxy 
will move back-and-forth between magnitude bins in the 
luminosity function. If the characteristic intensities, du- 
rations, or frequencies (duty cycles) of these episodes de- 
crease with redshift, the resulting effect will be to alter 
the shape of the luminosity function. 

In particular, if intrinsically low-luminosity galaxies 
are spending progressively more time in the state of ele- 
vated star formation (be it because the intense episodes 
are longer or occur more frequently) then the faint end 
of the LF will steepen; an increase in the intensities (i.e., 
star formation rates) of the starburst episodes will have 
a similar effect on the LF. At the same time, the bright 
end of the LF will remain constant if the characteristics 
of the starbursting episodes in the intrinsically luminous 
galaxies remain fixed with time. 

It is not obvious what mechanism could be responsible 
for the change in the duration or intensity of star forma- 



tion, but one possibility is that star-formation becomes 
more robust against self-disruption by feedback effects as 
their host dark matter halos accrete material with time. 
The lack of evolution at the bright end would then sug- 
gest that the mechanism depends on halo mass and has 
saturated for the more massive galaxies so that it cannot 
evolve further with time even though the dark matter 
halos themselves may still be growing. 

Scenario C: Evolution in the properties of dust? 

The third scenario we examine is linked to the possible 
evolution in the properties of interstellar dust in high-z 
galaxies. Even a small change in the amount of dust can 
have a strong effect on the observed UV luminosity of a 
galaxy while leaving its colors relatively unaffected. For 
example a decrease in extinction from E{B— y)=0.25 to 
E{B — y)=0.1 in a 2;~3.5 LBG would produce a six-fold 
increase in its M1700 — enough to match the evolution 
of the faint end of the LF we see between z~4 and z^2> 
— and yet would result in galaxy colors that still remain 
well within the LBG color-color selection criteria. Be- 
cause the luminosity change of the bright end of the LF 
is very strongly ruled out by the data, such evolution in 
dust properties would have to be differential with lumi- 
nosity, which suggests that such changes cannot be due 
to differences in sample selection between redshift bins, 
but could be due to real changes in dust properties. 

Again, it is not immediately clear what mechanism 
could result in a change in effective dust opacity in low- 
luminosity but not high-luminosity LBGs. If the dust 
evolution scenario is correct than it may reflect time- 
dependent changes in the properties of dust grains in the 
sub-L* LBGs or in the amount of obscured vs. unob- 
scured area visible in each LBG. 

While it is interesting to speculate about the nature 
of the underlying evolution of individiual galaxies that 
is reflected in the evolving LF, clearly, the luminosity 
function by itself is insufficient to discriminate between 
the possible mechanisms that are responsible. To under- 
stand what drives the changes we see, we will have to 
turn to follow-up studies that compare the properties of 
dust and star formation in high-z galaxies as a function 
of redshift and luminosity. 

6.4. The way forward 

A key result is that we have identifled luminosity and 
redshift as important variables in galaxy evolution at 
high redshift. We can use this fact to seek the nature 
of the underlying evolutionary mechanism by comparing 
diagnostics of dust, age, etc. as a function of L and z. 
While LBG follow-up studies to date have primarily fo- 
cused on luminous galaxies at z~3, now that we know 
that galaxy evolution depends on L and z, extending 
such studies as a function of luminosity and redshift pro- 
vides an attractive way to gain valuable insights into how 
galaxies form and evolve. 

As we have illustrated in S I6.3I evolution in the proper- 
ties of starbursting episodes or in the amount or distribu- 
tion of interstellar dust may underlie the evolution of the 
LF. One line of attack then is to compare the broadband 
spectral energy distributions of LBGs as a function of L 
and z. SED studies have already yielded insights into 
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the cxtincton, starburst ages and stcUar masses of LBGs 
(e.g., Sawicki & Yec 1998, Papovich ct aL 2001; Shaplcy 
et aL 2001) but by applying these SED-fitting techniques 
to different LBG sub-populations we can look for system- 
atic trends that may reflect the dominant evolutionary 
mechanisms. Evolution in the rest-frame UV-through- 
optical SEDs of luminous {L>L*) Lyman Break Galaxies 
from 2;~4 to z^3 suggests build-up of stellar mass and 
possibly a finely-tuned interplay between increasing dust 
content and star formation rates (Papovich et al. 2004) 
even in the observed absence of evolution in the lumi- 
nosity function; it will be interesting to see what SED 
analyses tell us about the evolution of sub-L* LBGs. 

Another approach will be to compare the detailed spec- 
tra of LBGs as a function of z and L. A composite spec- 
trum representing a 7?.~24.5 LBG has yielded detailed 
insights into the properties of LBG stellar populations, 
outflows, etc. (Shapley et al. 2003). Comparing com- 
posite spectra of LBGs of different luminosity and at 
different epochs may yield key insights into what makes 
LBGs different as a function of L and z. 

Yet another line of attack is to measure galaxy clus- 
tering as a function of both luminosity and redshift as 
this mcasiirement will let us relate the potentially time- 
varying UV luminosity to the more stable dark matter 
halo mass; while studies of the luminosity dependence of 
clustering have been attempted in the past (Giavalisco & 
Dickinson 2001; Ouchi et al. 2004b) they have either re- 
lied on vary small fields or on spectroscopically-untested 
selection techniques; the KDF is designed specifically 
with studying the dependence of clustering on LBG lumi- 
nosity in mind and we will attack this issue in KDF IV. 

There are two very important advantages that such fu- 
ture differential studies will have. Foremost, (1) from an 
experimenter's viewpoint, we now know that L and z are 
variables that affect how galaxies evolve. We can thus be 
confident that "varying" L and z will yield a "response" 
in galaxy properties linked to evolutionary mechanisms 
and that the lack of such "response" will equally impor- 
tantly rule out a candidate evolutionary mechanism. At 
the same time, (2) while such studies will likely build 
on previously-developed techniques as illustrated above, 
by comparing results as a function of L and z they will 
use these techniques in a an essentially differential sense 
thereby reducing our current reliance on theoretical mod- 
els or low-z analogs. Differential measurements are al- 
ways much easier and more robust than absolute ones, 
making such differential studies extremely attractive. 

We feel that important insights lie ahead using this 
differential approach and we will pursue such studies in 
the near future. 

7. SUMMARY AND CONCLUSIONS 

In this paper we have used our very deep KDF catalogs 
of UnGTZI color-color selected galaxies at high redhsift to 
construct the luminosity functions of UV-selected galax- 
ies at z^4, 3, 2.2, and 1.7. As we discuss in detail in 
KDF I, these catalogs use the very same color-color se- 
lection techniques as are used by Steidel et al. (1998, 
2003, 2004) to select their galaxy samples at these red- 
shifts. Moreover we use the same effective volume (V^//) 
approach to computing the galaxy luminosity function as 
used by Steidel et al. (1998) at z^3 and z~4. However, 
our KDF data select galaxies to TZiim=^7 — a magni- 



tude and a half deeper than that previous work - - and 
so allow us to probe the faint, sub-L*, end of the galaxy 
luminosity functions at these and lower redshifts. Our 
analysis probes the population to galaxies that have star 
formation rates of ^1 Mq/jt in the absence of interstel- 
lar dust and in our assumed (I^m, ^a, Ho) = 0.3, 0.7, 
70 km s~^ Mpc~^ cosmology. 

Spectroscopic redshifts for a large sample of galaxies 
to TZiim^27 would be observationally extremely expen- 
sive, and so the estimate of the luminosity function at 
these faint limits must at present rely on photometric 
redshifts or its cousin color-color selection. Several at- 
tempts to estimate the faint end of the LF at these 
redshifts have been made in the past. However, ours 
has the imporant and unique combination of using a 
color-color selection technique that is well-understood 
and well-tested spcctroscopically while at the same time 
drawing on a galaxy sample that's taken from a large area 
of 169arcmin^ thus giving good statistics and span- 
ning three spatially-independent fields — thus allowing 
us to control for cosmic variance. To our knowledge no 
other survey in existence has this combination of favor- 
able and important characteristics. We thus believe that 
ours is the most reliable estimate of the faint end of the 
luminosity function at these redshifts to date. 

We have carried out detailed studies to understand the 
potential systematic biases in our luminosity function 
analysis. We find that field-to-field variance or uncer- 
tainties due to fc-corrections do not significantly affect 
our results and find instead that the largest source of 
systematic uncertainty lies in the estimate of the effec- 
tive survey volume, 14//. We found that our results are 
robust to how we estimate K-// at z~3 and 2;~4; how- 
ever, the estimate of Ve/j introduces a source of system- 
atic uncertainty into the LF at 2~2.2 and z^l.7. Ad- 
ditionally, we also suspect that contamination by low-2; 
interlopers may be an additional source of uncertainty 
at z^l.7. Overall, we are highly confident of our LF es- 
timates at z~4 and 2:~3, we feel we can use the 2:~2.2 
results to place limits on the shape of the LF, and are 
not confident of the 2;~1.7 estimate. 

In light of the preceding discussion, the results of our 
analysis can be summarized as follows: 

1. The faint end slope of the LBG luminosity func- 
tion at Z'^3 and ^;~4 is shallower than the 
a=— 1.6±0.13 previously reported at z~3 by Stei- 
del et al. (1998) using Hubble Deep Field data. We 
find a^-lA3tom at z~3 and -1.26t[j;|g at z~4. 
While formally consistent with the Steidel et al. 
(1999) a, our more accurate, lower a may force 
a factor-of-two downward readjustment of many of 
the recent UV-based estimates of the density of star 
formation in the Universe at z^3 and above — an 
issue that we address in KDF III. 

2. We find strong evolution in the number density of 
faint (sub-L*) LBGs over the 600 Myr from z^4 
to z^3: there are 2.3 times more sub-L* LBGs at 
z~3 than at z^A. This result is statistically secure 
at the llfT level, and we believe it to be indepen- 
dent of systematic biases due to cosmic variance, 
sample selection differences, surface brightness se- 
lection differences, assumptions about the cosmo- 
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logical model, or differential fc-correction effects. 

3. While the faint end of the luminosity function 
evolves from z~4 to z~3, the bright end appears 
to remain unchanged. This differential, luminosity- 
dependent evolution is statistically significant at 
the 98.5% level, where the limitation in our confi- 
dence comes from the small number statistics of the 
bright end of the LF. An improvement in the level 
of confidence here will require analysis of Z'^A and 
2;~3 LBG surveys that are significantly larger than 
even the largest that have been studied to date. If 
the differential evolution is real then it may allow 
new, differential approaches to the study of Lyman 
Break Galaxies. 

4. It is not clear whether the evolution of the faint 
end of the luminosity function continues to lower 
redshift because of potential systematic biases at 
2:~2.2 and z^l.7. We find that our estimate of 
the z~2.2 LF depends on our assumptions about 
the make-up of the galaxy sample at this redshift. 
Despite these systematic uncertainties, we can nev- 
ertheless conclude that there are at least as many 
sub-L* galaxies at 2~2.2 as there are at z~3. 

5. At z~1.7, systematic effects make it diffcult to be 
confident of the reliability of our LF determination 
at that redshift. 

The two most intriguing results of the work presented 
here are the increase in the number density of sub-L* 
galaxies from z^A to z^3 and the possibility of differ- 
ential, luminosity-dependent evolution over that redshift 
interval. As we discussed in 14.31 and 15. II the increase 
in the number of low luminosity galaxies is a robust re- 
sult that is both statistically highly significant and at the 
same time unlikely to be an artefact of some systematic 
bias. The presence of differential evolution of the galaxy 
population is a less secure result (98.5% statistical prob- 
ability, § 15.2(1 that will need to be confirmed with much 
larger, shallow LBG surveys. However, we note that dif- 
ferential LF evolution is not unexpected given that there 
is no reason to think that galaxies across a range of UV 
luminosity — and so, to first order, a range of differ- 
ent star formation rates — are straightforwardly scaled 
analogs of each other. 

The evolution of the faint end of the population raises 
the intriguing question: what processes in individual 
galaxies underlie the observed evolution of the popula- 
tion? A wide range of possible evolutionary mechanisms 
may be at play, ranging from changes in the properties of 
star-bursting episodes that seem to occur in these galax- 
ies, to evolution in the amount or properties of interstel- 
lar dust. Discerning what mechanism is responsible will 
be important for our understanding of how high-redshift 
galaxies form and evolve. 

One avenue of attack on this problem is suggested by 
the fact that the evolution of the galaxy population may 
be differential with luminosity. If the population does 
evolve differentially with luminosity, then comparing the 
properties of Lyman break galaxies as a function of lu- 
minosity and (for faint LBGs) of redshift may point us 
towards the underlying mechanism. If — as some studies 



suggest (e.g., Giavalisco & Dickinson 2001) — UV lumi- 
nosity is a tracer of the dark matter halo mass, then we 
will first be able to link the differential evolution of the 
population to a mass scale. We are starting to pursue 
this line of attack with our KDF data (KDF IV). An- 
other set of insights will be possible from comparing the 
spectral energy distributions of LBGs or the details of 
their composite spectra; such studies have to date been 
focussed on relatively luminous LBGs at 2~3, where they 
have yielded insights into, e.g., the starbust ages, ex- 
tinction, and the state of their interstellar media (Saw- 
icki & Yee 1998; Papovich et al. 2001; Shapley et al. 
2001; 2003). The intriguing possibility of luminosity- 
dependent evolution of the LBG population opens up the 
possibility for such studies in a way that is differential 
and so largely independent of the systematics associated 
with using models or low-redshift analogs. We will ex- 
plore such approaches in upcoming work. 
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Fig. 1. — Colors of model galaxies in the UnCR, (left panel) and GTZI (right panel). The filled regions represent the color-color selection 
criteria used to select our samples of high-z galaxies (see KDF I and Steidel et al. 1999; 2003; 2004). The tracks show model colors of 
100 Myr-old (solid lines) and 10 Myr-old (broken lines) starbursts for three values of reddening each, E{B — V)=0, 0.15, and 0.3. Reddening 
generally increases from lower left to upper right. The points mark the location of z=1.7, 2.2, 3, and 4 on each of the tracks. 
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Fig. 2. — Comparison of our p{Tn,z) models with the observed redshift distributions of UV-selected high-x galaxies. The shaded 
histograms show the spectroscopic redshift distributions of Steidel et al. (1999, 2003, 2004), while the solid curves show our p(iri.z) for 
E(B — y)=0, 0.15, 0.3. The E{B — V')=0.15 fiducial model is marked with a heavier line and generally, in each panel the peak of the p(m, z) 
moves to lower redshifts with increasing E{B — V). These p(m,z) models arc for 7^=26.5 — a magnitude at which typical photometric 
errors in the KDF correspond to those in the observed spectrocopic samples whose photometry is shallower than ours. Here, results for 
our 09A field are shown but the other four fields give similar curves that show good agreement between our fiducial model model and the 
observed redshift distribution. 
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Fig. 3. — The fc-correction color term (top panels) and the distance moduli {bottom panels) for two representative galaxy SEDs. Locations 
of the z^A, zr^3, 2~2.2, and 2~1.7 sample bin centres are marked. The top panels show the color terms that are needed to transform 
observed G, GTZ, TZ, and / magnitudes to rest-frame 1700A. The top left panel shows the color term for a 100 Myr-old starburst with 
moderate dust attenuation, while the right panel is for an unobscured 10 Myr starburst. The errorbars show the FWHM redshift ranges 
(<52~0.3; see KDF I for details) spanned by the color-color selected samples, and the thick colored bands highlight the color term values 
corresponding to these redshift ranges. These plots illustrate that for the right choice of observed bandpass — namely / for z~4 LBGs, TZ 
for 2^3 ones, CR, for 2~2.2, and G for 2~1.7 — the fc-correction color is ^O. The bottom panels show the corresponding offsets between 
observed and absolute magnitudes (as defined in Eq. |B] — i.e., including the fc-correction color term), and illustrate that the redshift 
uncertainty for our photometrically-selected objects translates into only a small DM uncertainty, and hence into only a small uncertainty 
in the derived absolute UV magnitude of the object, A/1700 • 
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Fig. 4.— 



The rest UV luminosity functions at different redshifts. The solid black curves show the best-fitting Schechter functions, while 
the shaded bands show the corresponding 68.3% confidence region. As is described in the text, the errorbars include both the \fN statistics 
and a bootstrap estimate of field-to- field variations. The dashed fiducial curve simply reproduces the z^'i curve. The dotted curves show 
the GALEX rest-frame f500A LFs for comparison; the rightmost, fainter one is for zr^O and the leftmost, brighter one is for zr^l. As we 
discuss in the text, we regard the 2;~4 and 3 LFs shown here to be highly trustworthy, the 0~2.2 to represent a firm lower limit on the 
galaxy number density, but the z~1.7 to be questionable due to systematic biases. 
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The Schechter luminosity function parameters and their Icr confidence regions. 
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Fig. 6. — The cfTcct of changing cosmology on absolute magnitudes (top panel) and number densities (bottom panel) as a func- 
tion of redshift. The quantities plotted show the change incurred in transforming from the cosmology adopted in this paper — 
(f^M I ^^Ai ^^o)=(0.3,0.7,70) — to the two alternative cosmologies labeled in the plots. The values of M and (p in the alternate cosmologies 
are denoted with primed quantities while in our cosmology they are unprimed. While the change in M or (p at any given redshift can be 
quite large, the relative change from redshift to redshift is small, ensuring that — for reasonable cosmologies such as those considered here 
— any evolutionary trends seen in the LF are not subject to the assumed cosmology. 
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Fig. 7. — Dependence of the LF on the Ve// model. Variation in the derived LF for the range of considered Ve// models. The 
points and black curves show the results of recomputing the LFs using the seven different models of Vef / that result from seven different 
E(B — y)values. The gray curves are the z^3 results replotted in the other redshift panels. See text for more details. 
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Fig. 8. — Dependence of the galaxy number density at Mi7oo=— 20 on the V^ff model. The black symbols show 100 Myr-old starburst 
models, and the gray ones are for 10 Myr-old models. 
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Fig. 9. — Comparisons with other surveys. Rest-frame UV luminosity functions measured in the FORS Deep Field (FDF; Gabasch et 
al. 2004) at and 2~4 and in the Subaru Deep Field (SDF; Ouchi et al. 2004a) at 2~4 are shown. Also shown are the results of the 
Steidel et al. (1999) analysis of the HDF. The Lyman Break Galaxy LF is shown as filled circles with errorbars that include an estimate 
of the field-to-field variance determined via bootstrap resampling. The FDF measurement is based on a single small field and does not 
include an estimate of cosmic variance while the SDF measurement is based on a single very large field and so is unlikely to be strongly 
affected by cosmic variance. The LBG luminosity function is in very good agreement with that found in the smaller-area FDF at both z^S 
and 2;~4. However, while the KDF and the FDF agree with each other, they both disagree with the SDF measurement at 2:~4. 
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Fig. 10. — Comparison of the z~4 and z~3 luminosity functions. The top panel shows the z~4, LFs together, with the z~3 Schechter 
fit marked with a heavier line. The bottom panel shows the quantity which measures the fractional deviation of the data (at z~4 or 3) 
from the z~3 Schechter fit: perfect agreement between data and the z~3 fit would put the points on the horizontal $=1 line. 
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Fig. 11. — Comparison of the 2~4, 3, 2.2, and 1.7 luminosity functions. As in Fig. llll but adding the 2~2.2 and 1.7 LPs. As is discussed 
in the text, the z~4 and z-^S LPs are robust to the details of the analysis procedure, whereas the 2~2.2 and z~1.7 are less so and are 
subject to potentially strong systematic uncertainties. 
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Fig. 12. — The luminosity-dependent evolution of the number density of galaxies. The horizontal axis shows the amount_of change in 
the density of luminous galaxies (Mi7oo<— 21) parametrized as the LF as the ratio between the average normalized densities <E> at zr^A and 
z^S. The vertical axis shows the same quantity but for low-luminosity galaxies, Mi 700 >— 21. The straight lines represent three fiducial 
cases: no faint-end evolution (horizontal line), no bright-end evolution (vertical line) and equal evolution at the bright and faint ends 
(diagonal line). The point and ellipses represent the amount of evolution from z~4 to 2~3 and the associated 68% and 95% confidence 
regions. There is substantial evolution of the low-luminosity population from to 2~3 that is statistically significant at the 98.7% level. 
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Table 1. Parameters of Schechter Function Fits 



z 


Steidel type^ 


^-'^1700 


log 4>*^ 


a 


1.7 
2.2 
3 
4 


BM 
BX 

C, D, M, and MD 


-19-80lo.26 
-20.60t°-j^ 

-20.90+0^4 

-21 00+°-*" 

Z1.UU_Q 4g 


-1 77+0.11 
'-0.11 
r, c;o+0.20 
^■'-''^-0.26 
r, 77+0.13 
^- ' ' -0.09 
_o Q7+0.21 
-0.33 


-o.8it°:- 

_1 43+0.17 
-1 26+°''° 

^•^"-0.36 



''In the nomenclature of Steidel et al. 2003, 2004 
''In units of Mpc~^ 



Table 2. Number density log 0*^ in different fields of the survey 



z 


combined fields'" 


02A 


03A 


03B 


09A 


09B 


1.7 
2.2 
3 
4 


-1 77+0.11 
' '-0.11 

r, t.r, + 0.20 

^•''^-0.26 
r, 77+0.13 
^- ' ' -0.09 

'J-'-" -0.33 


-1.92 
-2.59 
-2.79 
-2.97 


-1.80 
-2.41 
-2.69 
-3.24 


-1.80 
-2.44 
-2.88 
-3.27 


-1.70 
-2.65 
-2.77 
-3.16 


-1.74 
-2.55 
-2.74 
-2.94 



^In units of Mpc ^ 

''Prom fit to the full data as reported in Table [Tl 



